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ABSTRACT 


The present thesis will allow Users to analyze rectangular 
plates with or without a rectangular hole with its edges parallel 
to the edges of the plate, by using an approximate grid analogy. 

The method is completely computerized and the user need 
only specify properties, geometry, loading and kind of output de- 
sirable. The computer programs then will generate a deck of cards, 
that form a complete STRUDL input. STRUDL itself is an equivalent 
of STRESS but with mich more useful capabilities. 

Initially it was intended to make everything an one step 
process, i.e. to read input in a free-field form, generate all 
necessary statements of STRUDL, store them on a file and let 
STRUDL read the file and give the desired output. liowever due to 
several delays in the Civil Engineering system Laboratory, STRUDL 
is not yet available and CDL (Command definition language) is now 
in the process of experimentation. In consequence the user has to 
supply his input according to specific format and the linking be- 
tween this subsystem and STRUDL is not automatic. In other words 
the process is a two-step one and the user has to take the output 
of this program and use it as a STRUDL job. 

The solution to the flexural plate problem contained herein, 
is extremely general regarding boundary conditions and loadings 
so that almost every engineering preblem of this nature can be 


taken care of, 








The particular connecticn of the present job with Naval 
Architecture is that selutions to bottom and side plating with 
hydrostatic loads, bulkheads with openings and decks with cargo 
holes will become readily available with ea minimum of human effort 


and intervention, 
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CHAPTER I 





The develoument of the large digital computers was a decie+ 
sive step towards the solution of many engineering problems. Indeed 
during the past few years people looked back and reoxar 


but it was impossible to carry it out by hand. it was a sau. 
simple matter to solve all these problems then by a digital com- 
puter. 

Among all these other problems one is the structural analy- 
sis of plates a special case of which is the rectangular plate in 
flexure. Several scientists and engineers have done a marvellous 
job in solving special cases of it. 

it seems that there would be two distinct approaches in 
coding the solution of the problem into a commuter: 

(4) Take each special case and code its exact solution. 

(44) Invent a unified approach applicable in all or most 
special cases, 

The first approach requires a treserndous ammount of work. 
Besides it could never be called general. For, no matter how many 
special cases one solves,there will always be some other special 
case left out. 

The second approach seenus quite attractive, with the only 
disadvantage that it is an approximate one, However, in cases where 
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an exact solution is not known, an approximate one is best than 
nothing. Eesides, very often engineers are interested in getting 
a first estimate of what is happening in a system. 

in the structural analysis of rectangular plates, the grid 
analogy which will be deseribed in dotail in the next chapter, is 
among tke best approximate methods to attack the general problem. 
in addition since several methods to solve skeletal structures 
have now been developped, one has only to transform the plate 
into such a structure, that existing computer forsulations can 
golve. 

The geal set forward by the author has been to develop 
a set of computer programs, able to reduce any rectangular plate 
to its skeletal equivalent and then solve it as such. Although 
the sample problem included in the present work is a very simple 
one, the capabilities of the prosrams develcped are quite general. 
In fact the choice of the sample problem, made strictly because 
{ts solution already exists, is a very unfortunate one in illustra-~ 
ting the application of the method into the more general problen. 

The case of a rectangular hole on the plate was included 
here, due to the interest of Naval Architects on the effect of 
cargo openings on the strength of ship's decks. Also maneholes 
on vertical bulkheads forming bounderies of fuelleoil tanks con- 
stitute another field cf application. 

What really makes the method powerful is its generality 


regarding the boundary conditions and the loading. 


indeed, roint supports at any point over the area of the 
plate, line boundary conditions over all or any part of the edges, 
concentrated point loads, uniferm or linearly distributed line 
and area loads or any combination of them, are sll cases tiat can 
be taken care by the programs with a minimum of effort. in addition 
the ability to get any kind of output desired might be very useful 
in some Cases. 

In the following chapters the analogy will be described 
and the method of its computer application will be explained. in 


addition directions will be given to the users on how to use the 


programs. 
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CHArTER IT 


2ele= General 

Perhaps one could say with complete justification that the gene~ 
yal problem of a plate in flemre, has been solved, at least in principle, 
indeed the differential equation deseribing the phenomenon is known for 
many years ard it ia only necessary to cmnsider the loading an! bournlary 
comitions applicable in sach particular problem, However this last. step, 
namely, the solution of the differential equation for a given set of 
boundary conditions and a given loading, senetines presents inmamerable 
difficulties from the analytical point of view. Bven in the case where 
a solution can be found in a reasonable mmber of steps, it is very 
often expressed in the forn of an infinite series, which especially 
in the case of moments and shears converge rather slowly, a fact thet 
reduces their value. 

It is the above difficulty that sade engineers search for cther, 
perhaps less accurate but more easily developed solutions of the plate 
problem. There have been three approaches that yieldapproximate solu» 
tions to the plate stresses problen. 

(4) Finite differences methods, 

(14) Grid analogies. 

(414) Finite element methods. 

Cf these, the finite differences approech, depending uron the 


Se 


boumlary conditions, very often requires fictitious extensions of 
menbers or consideration of members with fictitious properties. This ma- 
kes a rather complicated task the coding of the general rectangular plate 
solution to a digital computer. There is also some suspicion that this 
methed gives results deviating from the exact ones (those that they 
exist) more than those computed by the other two methods. 

The finite element method seems to be quite promising not only for 
the rectangular plate case but for shells of arbitrary shape. Conside- 
rable mmber of people is working on that mothod, in the clvil Engineering 
department of H.I.T., but there will be some time before this methed bee 
comes available to users. 

The grid analogy was selected by the author for the following 
reasons. 

(4) It 4s rather simple and gives very satisfactory results 
regarding moment intensities and deflections as it was indicated in (4), 
(2) and (5). 

(44) It can be rather easily programmed to a digital computer. 

(444) The solution to the equivalent grid can be very easily 
found by the existing STRESS program or by its new equivalent STRUDL, 

(7) ,€8) (9). 

(iv) The present work intended primarily for the field of Naval 
Architecture as well as other professions interested in structural 
engineering deals almost exclusively with rectangular flat plates or 
plates with small curvature which can be neglected, without considerable 


reduction in accuracy. 
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2.2. Description of grid analogy 
The grid analog of a flat plate has been represefitei differently 


by different people. Lightfoot in (3) uses a square gridwork consisting 
of grid beams with flexural rigidity and a set cf torsion bars with tor- 
siomal rigidity with no diagonals. Efimba in (6) in a program written 
to analyse bridge decks only, uses Tamberg®s method of analysing bean 
and slab type bridges as torsionally weak grids (10). 

The appreach taken by Hrennikoff in (2) and Yettram and Husain in 
(1) 4s, that the plate can be represented by a plane structure composed 
of articulated bars, jcined by pins normal to the plane of the structure, 
The pattern has to repeat itself throughout the area of the plate and 
the outline and external dimensions of the plate should be the same 
with the corresponding ones cf the equivalent structure, Repetition of 
the pattern does not however mean that the size of the unit has to ree 
main the same. In fact by changing the size of the unit we can satisfy 
ary existingdeviations from the simplest case of a sinple rectangular 
plate, like differences in thickness, holes or differences in material 
properties, 

Hrennikoeff set up the problem in a very detailed fom andi 
developed the equations that will define the properties of the elenents 
of the unit. The restrictions that he imposes in crder to solve for 
the properties of the members are, that the frame work has to deflect 
properly: 

(4) When subjected to pure bending in the X direction with 
no bending deformations in the Y direction. 

(44) When in pure bending in the Y direction with no deforma- 


tion in the X direction. 
(iii) When subjected to pure torsion on the X and Y planes. 


These are general restrictions applicable in any case and are 
useful in examining the suitability of any particular unit to be 
investigated. The third condition falls out whenever the pattern has 
the same form in the direction of both axes. 

Subsequently Hrenikoff investigated two patterns as shown in 
fig's (2.1) and (2.2) and proceeded to develop their properties. 

Of these two patterns the first one restricts the solution 
to values of Poisson's ratio of 1/3 only. 

The second one although valid for 
any Poisson's ratio has the disadvan- 
tage that is a square one. This causes 


a serious problem in the case that 


A 
we cannot satisfy the boundary don- 
ditions with a square pattern. We 
ae 
FIG. 4 also notice that for each square we 
SIMPLE SQUARE PATIEQAN have 10 slender members, which would 





make a hand-solution more difficult 
or a computer solution lonzer. 

A rectangular pattern with diagonals 
only is much more general and simple. 
This type of pattern has been very 
recently investigated by Yettram 


and Husain in (1). Their pattern 


Det 


FIG. 2.2@ SQUARE PATTERN 
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Looks to be the best of the existing ones because it is simple, easily 
developed and coneiderably general. Due to the fact that the seth 


is relatively new andl espesially useful in other kinds of problems 
(those satisfying the Langrangian equation like thermal stresses in 
& plate of special cases of fluid flow), it is believed that a rela 
tively detailed discussion ot it mst be given. This will be dame in 
the next chapter. 
fo3. Sectangular pattern crid anglogy 

The id@a of the model is that we can substitude a rectangular 
network consisting of slender sembers only, in place of 4 particular 
plate clement. Fig. (2.3a) shows a rectangular element of plating and 
Fig.(2.3b) its equivalent grid element. Tne unknown properties of the 
grad element members are the flexural rigidity of the side beans I, 
and I, , the torsional rigidity of the side deans ee and ee A 
and the flexural rigidity of the diagonals I, . That. gives 5 unknewn 


G 
properties i.e. I, , I, + la, S and oie 


The eyjuations giving the above 5 properties can be developed 





by equating the rotations at the nodes of the grid element with those 
of the plating element of equal size. Although the methed doesn't zive 
the crogs-sectional area of the beams, it gives all the necessary ine 
formation to make the solution of the grid frameework possible. Of 
course the plate and the grid elements must be subject to statically 
ejuivalent torgues and moments, 

Considering first the plating elenent we examine three cases. 


(1) “qual and oppesite bending moment intensities on the x sides 





(c) a (d) 





FIG. 2.3 PLATE AND GRID ELEMENTS 
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(41) Equal ard opposite intensilies on.y sides, 
(411) Self-eguilibrating torque intensities along 4 edzes. 


Fig's (3d), (3e) and (3g) show clearly the toment and toraue 
distributions and the corresponding angles. The values of the angles 
in terms of the moment intensities can be found by differentiating the 
expressions for the deflections, (11). This was done in (1) where 


they came up with the following expressions for the anzles, 
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For the grid slement loaded as shown in Fig's (34), (3f) and 
(3h), Ref (1) gives the flexibility matrix (shown here in Fig 2. 4) 
and concludes that only 2 out of the le equations are indevenient when 
we have only side bending moments or self-equilibrating torques. 
By selecting properly these 2 equations in each case we come up with 


the following expressions for the grid angles. 
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By equating the rotations for the plate and the grid clements 
we can get the following equations giving the required properties of 


the grid member. [The method iS explained in more detail in (1)], 





teem: Se, m (2,12) 
pees <a (2.13) 
Ty, Kod, (2.14) 
Gig KGa 
a Cees ae (2.15) 
ice. (apo Ee 





ee, ie pai) 








14, 


Where k is the ratio of the side beans 
L ig tne Poisson's ratio 
h is the thicknesa of plating 
& is tae Soung's xodulus of elasticity 
& le the modulus of rigidity 
® ise the diagonal te side beam ratio 
For future use we will transform the abcve relations by using: 


a t= = (2.17) 
g(1+p) 
into the following fora: 
Qsigha ae (2,18) 
ge = S \ 


T.= FAGH KE) (2,19) 
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3. (2.20) 
Se ee ee (2.22) 
a al ad ule (2,22) 
where - Serene (2.23) 


is a constant for a ziven isotropic plate of constant thickness. 


The potentialities of the above equations ere inuediately 
apparent. They can be used for the solution of any plating problem, 
with arg boundary conditions ami load and any mmber of beles provided 
that we can form a grid fine enough so that it will represent a parti- 
qilar plating element, i.e. provided that we will satisfy equations 
(2.18) through (2,23) at each particular element in the division of 
the plate. Also by forming a proper grid we can take care cf non- 


isetrepic plates. 
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In the following chapter some limitations and principles of 
appiication of the above nethed will bo discussed before the computer 


approach is explained. 


24, mitations and oringiples of apulicaticon of Grid method 

Expressions (2.18) through (2.23) might give, under special 
clrounstances negative values for the I's or J's, Altscugh Ref. (1) 
clains that this is not a probles, existing computer programs fer the 
solution of a grid expect. positive moments of inertia, This linits 
the range of applicability of the method in some mamner, 

The range of the Foissen’s ratio fer engineering pur;oses 
varries between .15 and .3 . Looking carefully into equations 
(2,18) through (2.13) we conclude that in order to assure positive 
values of all J's and I's the value of K wust be: 


Vr <a VR (2,24) 


ee ee eee 2 (2425) 


For higher values of . than 1/3, negative properties are 


unavoidable. fewever the authors of Ref,(1) state that for good roe 
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sults the value of k mst be between 2 ard 0.5. 

The computer programs will be such that will automatically 
form a grid with k's (if more than one) as mich close to unity as 
possible. It is obvious here that the best results are obtained for 
k=] but sometimes the boundaries are such that 4t is quite difficult 
to have a square pattern. 

In the actual application of the method the following points 
mist be taken into account . 

(i) The properties found above give the contributions of only 
one element in the division of the plate. Thus for the side beam we 
have also to add the prorerties of the neighbouring clewents. Of 
course this is not tre for the diagonals and the members lying along 
the eiges of the plating or along the edges of any existing holes. 

(ii) The boundary conditions along the edges of the plating 
will be considered as applicable on the nodes of the grid at that 
neighbourhood, 

(441) Any isolated point supports in the interior of the plate 
are considered applicable on the nearest node, 

(d) The loading is considered as acting on the nodes or on 
the members of the grid. 

The exact manner in which the loadings and the boundary condi- 
tions are treated will be explained in the description of the parti- 
cular SUBROUTINES that handle them. Of course most of the assumptions 
made here require that we have a fine enough grid so that the effect 


on changing the problem from the criginal one to another problem will 
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CHAPTER ITI 


A DRIEF DISCUSSION OF ICES 


3.1. General 

Although this chapter seems fairly obvious to people in the 
Civil Engineering Department of M.I.T., related to computer work, 
it is strongly believed that at least a brief description cof ICES mst 
be siven here for people not familiar with it. 

ICES (Integrated Civil Engineering system) is a computer systen 
designed for the efficient solution of a variety of Civil Engineering 
problems among which one is the structural Analysis Problem. 

Although it looks as the system is applicable only to Civil 
Engineering problems the generality of some SUSSYSTEMS of it makes 
them quite useful to any other field of Engineering, 

The system has been espesially designed in such a way that the 
user (A user is a person that uses programs written by other people 
called the programmers, for engineering computational work) communi-~- 
cates with the computer by informing the computer on the nature cf 
the problem and recommending the steps to be taken for the solution 
of the problem, In additien he doesn't have to previde the input at 
ary particular fcrmat. In that way the engineer instructs the computer 
instead of following instructions imposed by the machine. 

ICES has been made to fit into the 360 ILM machine and uses 


as a programming lancueage, ICEIRAN, ICEIRAN (ICES, FORTRAN) is a com- 
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bination of FORTRAN IV slus vowerful new capabilities useful in 
performing engineering work, All ICES programs are suprosed to be 
written in ICHSRAN. For the murposes of this work it suffices to say 
that a program in ICETRAN is crecompiled by a special precompiler 
(the ICETRAN precompiler) into a FORTRAN program which is finally 
processed by the FCRTRAN compiler. 

A SUESYSTSH in ICES is a set of pregrams that can be used 
independently or ccllectively to solve engineering sreblems. One of 
the subsystews is STRUDL (STRUCTURAL ANALYSIS & DESIGN) which includes 
all the capabilities of STRSSS plus new ones. In particular the work 
done herein will hopefully initiate another subsystem ealled PLATES, 


Which eventually will be able to solve al]1 kinds of problems related 


to plates. 





There have been several capabilities of 10S and ICETRAN that 
simplify very much this work and in addition make it quite general. 
These capabilites are: 

(a) 
problem has been to reduce it into en equivalent grid and then let 


ms. Since the apvroach to this 





STRUDL solve it, the interaction between subsystems is very important. 
At the time that this work was untertaken this feature was not avai- 
lable and that is why it was decided to make the cutmut of this 
program a deck of cards that form a complete STRUDL input. However 


it is a simple matter to modify this program to assure one-step 
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processing when interaction between subsystems becomes available. 
(eo) The ICES executive program. It is the ICES executive system 


that gives the programmer the capability of allowing the user to write 
his input in an English-like fora. The executive process consists of 
reading the input, finding the command in a dictionary, storing the 
data in the proper place and executing associated subroutines, 

The rules for defining commands and associated crocedures are 
handled by the command Definition Language (CDL) compiler. More ine 
formation about how to define commands, what the ICEX(ICES SXE&CUTIVE) 
input structure looks like ete., may be found in (9), chapter 3, 

(c) Dynamic memory allocation 

The characteristics of the equivalent Grid might change from 
one proble= to another. We cannot then say that we will arbitrarily 
set the number of nedes to 500 for examble, since this will depend on 
the kind of the problem in general and the accuracy that the user 
desires. That. indicates that the arrays to be used throughout the 
probles: are not of a given size but their size varries from problem 
to problem. As it is well-known in FORTRAN one would have to indicate 
the size of the different arrays at tne beginning of the progran. 
Then if the problem is small (relative to the defined size of arrays) 
menory storage is waisted. If the problem is large it cannot be exe- 
cuted, 

By using the dynazic mexory allocation in ICEIRAN the programmer 
ean define the size of the arrays during the execution time. This in 
tarn has the effect that at no time the menory is waisted, 

In addition to being able to define the size of the arrays 
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during the execution period ancther advantarce of the dynamic memory 
allocation is thse fact that secondary memory is considered a logical 
extension of the primary memory. Thus at ary tine that the primary 
monory is filled up, the secondary memory automatically is used to 
store additional information. 

the important thing then for the user is that the size of 
his program is not restricted by the primary menory size. what makes. 
the dynamic memory allocation powerful and nice is the fact that 
neither the programmer nor the user has to worry about size of arrays, 
provided that thay are dynamic and their size has been defined by 
an arithmetic expression. 

(d) 

The last extremely useful feature,of ICETRAN, regarding this 





work ig the organization of secondary menory, Since the data generated 
by the PLATE program will be finally used by STRUDL it looks quite 
right to write all of it on a file and then let STRUDL read that 
file and proceed to the solution. 

This has not been done at this stage because of the following 
reasons. 

(i) When this work was done STRUDI, didn't have the capsbility 
to accept infermation from a file in an intornal form. 

(14) The system ROUTINES that can be used to write information 
on a file cannot do it under a specified FORTRAN type FORMAT, a fact 
that would allow the writing of information on a FILE in an external 


forn, 
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(444) There are statements in FONTRAM that ean be used to write 
information on a PILE under a szecified FORT statment but the sige 
of the file must be specified in advance (Like the DIZNSIGN statenent) 
which makes them inadequate for our case, 

However when the capabilities (a) and (b) will become aimilable 
3t will be a simple matter to modify the programs in such a way that 
information will be directly written on a file wrich STRUDL will be 
able to read. 

This eoneludes the brief discussion of ICES given in this 
chapter, The purpose of it, was to give to the none-familler reacer 
a slight idea, why this work has been done in ICETRAX, A complete 


description of ICKS is given in (9), 
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CHAPTER IV 


CUMFUTER PROGRAMS 


4.1. Introduction and capabilities 


Eefore a general description of the computer programs is given, 
it is necessary to give arather detailed description of the capabi- 


lities of the new subsystem. 


b 


0 b 


F\ Lo = LOCATION AND ORIENTATION OF AXES 


Fig. (4.1) indicates the orientation and location of the axes 
of coordinates relative to the plate. The 2 - axis is positive upwards. 
This subsystem will be able to handle constant thickness rectangular 
plates of canstant Young's modulus and Poisson's ratio. The capabilities 
of the subsystem are: 

(i) Number of holes: There can be either no hole at all or a 
rectangular hole with its edges parallel to the edges of the plate, lo- 
cated anywhere on the surface of the plate. 

(ii) Boundary conditions: The boundary conditions that this work 


can handle can be divided into two categories. The first category 


Hi iH] 
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includes line boundary conditions (hinged free or fixed) arround 311 

or any gortion of the edges cf the plating, Needless to say that any 

one of the four edges cf the plating might be hinged along one vortion 
of it, Sree along another and fixed over the rest of it, The second cate- 
gory includes simple point suerorts over the entire area of the 

plate except the hole if there is one. No line supports are allowed at 
any other locations except tre four edges. however any line surrort, 

in a position other than the edges, can be perfectly well treated as 

a set of equally spaced point supports. 

(144) Loading: The loading can be of one of three types. A 
vertical force along the z - axis or a moment the vector of which 
is parallel to the x or y axis. If the vector of the moment has an 
orientation other than parallel to the x or y axis it can be treated 
by considering its components. In any of the above three cases the loa 
ding can be a concentrated one acting on a point, a uniform one along 
@ straight line or a lineer one along a straight line. Forces not 
parallel to the z- axis or bending moments having a component parallel 
te the 2 - axis connot be treated since this requires a different 
grid analogy. 

(iv) Gutput: Moments forces or reactions at specified points, 
along straight lines, or over the entire area of the plating can be 
requested, 

A more detailed description explaining how each case has been 
taken care of, will be given in the documentation chapter of the vorogramns 


that perform szsecific tasks. 
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4&2. A macroscon the pram 

The entire work can be legically divided into 4 steps, namely: 

(1) Formation of the equivalent grid . 

(42) Reduction of the boundary conditions inte ejuivalent 
grid supports. 

(141) Reduction of the loads into equivalent grid loads acting 
on the joints or the members of the grid. 

(iv) Reduction of the specified output into an equivalent 
legitimate grid output. 

Each of the 4 steps outlined above is performed by a set of 
subroutines. Some of tho subroutines are equally applicable in both 
eases Where a hole does or doesnot exist. Whenever this is not 
possible a different subroutine is used for each case. The USER 
however doesnot have to worry about these details. If he specifies 
that there or there is not a hole the proper subroutines are exe- 
cuted automatically. 

In the next chapters a macroscopic description of each of 


the above 4 steps will be given. 





The first step in the solution of the problen is the gone~- 
ration of the equivalent grid. This step is a combination of zeo- 
metry and particular boundary cenditions in order to produce the 
proper JOINT COORDINATES, JOINT RELEASES and MEMBER INCIDENCES 
commands in STRUDL. This part of the job is handled by the following 
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subroutines: 
DIFOHO (Division Points ~ Hole) 
DIFONH (Division Points ~- No Hole) 
GRID (Geometric Formation of Grid ~ no hole) 
MUGRID (Geometric Formation of Milti Grid - there is a hole) 
MEINNH (Member Incidences - Ho hole) 
HEINHO (Member Incidences ~ Hole) 
MEMPRO (Member properties) 
JOLNRE (Joint Releases) 
A generalized flow-diazram is given in Figy.51. explaining 


the sequence and interaction of this phase as well as the other 


three phases. 





The boundary conditions of the plate are transformed into 
an equivalent set. of supperts for the grid by the same subroutines 
no matter if there is or there is not a hole, 

The specific routines,that do the job here are: 

POBOC (Point Boundary Conditions) 

LIBCC (Line Boundary Conditions) 


4.5. Loading 


There are four subroutines that perform the transformation 


of the plate loading into an equivalent grid leading. These are: 


LINELC (Line Loadings) 








FOLNLO (Point Loadings) 
LINE (A logieal extension of LINELO to get arround the difficulty 
of the :RECONPILER to treat large mumber of syxbols) 


AREALO (Ayes Loading) 


Ne6. Cutout 
Theres is enly one subroutine (OUT) which in conjuction with 
POIRLO and LINELO is taking care of the output. 
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CHAPTER V 


PROGRAM DOCUMENTATION 
Sel. Seneral 

This chapter contains adequate information on how the dif- 
ferent programs function and how they are linked to each other, 
This information is — especially for other programmers 
intending to modify a particular subroutine, In fact this chapter 
is mostly useful to other subprogrammors but it can also be of 
some help to the users of the program. 

Especially for users the chapter on @ senple problem and 
form of data is extremely useful, in guiding them how to write 
their own input. 

In writing the documentation of the SUBNCUTINES the author 
conformed with the recommended methods of (9). Particular attention 
hes been paid on including all the fine points of the prograns 
into either their long description or their flow diagrams, However 
if any point is not made clear in the documentation cf a program 
the reader is referred to Appandix B which contains all the ICEIRAN 
listings. 

Some ceneral points applicable to the entire documentation 
program are: 

(4) All COMMON*'S are the same throughout all SUBRUTING. 

(41) All global dynamic ARRAYS are defined in eitaor DLPONH 


or DIPCHO and are destroyed by the MAIN FROCRAM, 
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(444) All LOCAL DYNAMIC ARRAYS are defined in the particu- 
lar SUBRCUTINE that uses them and destroyed before a return is 
executed from the subroutine. 

(iv) All transmitted dyn. arrays are defined in the CALLING 
program and destroyed before a return is executed from the calling 
preorram. 

(v) Sinee it 4s impossible at this time to have 1/0 state- 
ments with dynamic array elements, all cutrput is written by first 
storing the dynamic element into a nonedynamic one and then using 
I/O statements with this last clenent, 

In the following chapters a description for each subrou- 
tine is given. Appendix A contsins the flow diagrams, Appendix 5 
the Icetran listings and Appendix E a dictionary of all names in 


COMMON, 








Bhs 


Se2 =SUBROUTIES DCIPORM 





1 - NAME « Subroutine DIPom 

2 » AUTHOR = John Kolliniatis 

3 = DATE « May 7, 1966 

& » TYPE ~ Specifie 

5 @ SORT DISCRIFTION - This subroutine is used in the case that 
there is ne hole, in order to compute the mmber of division 
points along each side for a srecified maximun of Joints. 

6 = LANGUAGE ~ Ieetran 

7 ~- FLOW CHART ~- As shown in Fig (45.2) Appendix A. 

8 » USAGE ~ CALL DIPONH with the necessary input values in COMMON, 
Input variables are: 
VELEN = VELENI = Span of plate 
HOLEN = iOLEM s Width of plate 
MANUJO = Desired maxirun rmumber of joints, 

9 » OUTPUT ~ Cutput in COON includes: 
NUHOJ1 = Number of division points along the x-axis 
NOVEJ1 = Number of division points along the y-axis 


NJ 

NUKOMB 

HOVEMB With the same meaning 
HURIDI as in subroutine DIPOHO 
MILEDI 


AGe LENCGTin 
Sources Program - 60 Ieetran instructions 
Object Program - 67 Binary cards 

die CALLS - None 

iZe AFPLICATIONS ~ SUBSYSTEM PLATES 

43 + LONG DISCRIPTION + This subroutine computes the division pts 
along either side of the plate in such a way that the spacing 
of them will be equal. Since the values found will be fractioe- 
mal it truneates the values computed to the lower integral. 
Subsequently it examines whether it is possible to increase 
the division points aleng only one side so that we approach 
more to the specified number of division pts. 
in addition, since this subroutine is the first cne to be 
executed if there is no hole it. defines all the global dyna- 


mic arrays. 
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SUBROUTINE DIPOXO 

NAME = Subroutine DIPOHO 

AUTHOR - John Kolliniatis 

DATE - April 20, 1966 

TYPE - Specific 

SHORT DESCRIPTION - This subroutine is used in the case of a 

rectangular plate with a hole and a specified maximum number 

of joints for the equivalent grid, in order to find the sub- 

division points for each of the eight sub-panels that the 

plate is considered to be consisting of. 

LANGUAGE - Icetran 

FLOW CHART - As shown in Fig. (A5.3) Appendix A. 

USAGE ~ Call DIPOHO, with the necessary input values in COMMON, 

Input variables are: 

HOLEN = Length of plate along x-axis 

VELEN = Width of plate along y-axis 

HOLENI = Distance along the x-axis between the center of co- 
ordinates and the lower left corner of the hole. 

VELENI = Distance along the y-axis between the center of 
coordinates and the lower left corner of the hole 

HOLEN2 = Length of hole along the x-axis. 

VELEN2 = Width of hole along the y-axis 

MANUJO = Desired maximum mimber of joints for the equivalent 


grid 
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9 —- OUTPUT - Cutput in COMMON includes: 

NUHOJ1 : Number of division points along the x-axis of sub- 
panels A, Dand F (seo Fig. 5.1 and 5.2). 

NUHOJ2 : Number of division points along the x-axis of sub- 
panels B and G. 

NUECQJ3 : Number of division points along the x-asis of sub- 
panels C, E and i. 

MUVEJi : Number of division points along the y~axis of sub- 
panels A,B and C. 


NUVEJ2 : Number of division points along the y-axis of sub- 


panels D and &. 
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Number of division points along the y~axis of sub- 


panels F, G and H. 


: Number of joints. 


Number of horizontal members 
Number of vertical members 
Number of right - diagonals 


Number of left - diagonals 


; Total number of members 


MUECGJ1 = 4 
NUHOJ2 = 4 
KUHOS3 = 2 
NUVEJ1 = 4 
NUVEJ2 = 3 
NUVEJ3 = 3 
NUHOMB = 53 


NUVEMB = 52 


NURIDI = 43 
KULZDI = 43 
Wie te 
MEMBER 
4 NJ = 62 
ot Cex Rigur 
Per Gee ee DIAGONAL 
FiG-S5.Q- SCETCH IeLYUSTRATING NOMECLATYRE 
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10 = LuNGiz 
Source Program - 70 ICETRAN instructions 
Object Program + 50 Binary cards 

il. CALLS - Hone 

Ze APPLICATIONS ~ SUBSYSTEM PLATES 

430 LONG DSSCRIFTION - {[nit4ally this subroutine exanines whether 
it is possible to select values for the division points of the 
sub-panels in such a way that any crid element is a square one. 
Since these results are in fractional form in general a trunca- 
tion is made to the lowest integral value. In oase the hole is 
very close to one of the sides of the plate, the program ovor~ 
rides the restriction imposed by the USER by his specification 
of MANUJO, and arranges the division points in such a way that 
ensures consideration of the hole effect. Of course in this 
ease WJ is larger than MANUJC, 
In addition since this subroutine is the first one to be ex- 
ecuted in the case of a plate with no holes it is used to de« 
fine all global dynamite arrays. 


& 4, SUBROUTINE GRID 

4 = HAMS ~ Subroutine GRID 

2 = AUTHOR = John Kolliniatis 

3.~ DATE » May 7, 1966 

& » TYPE = Specifie 

§ » SHORT DESCRIPTION - Tnis subroutine is exeorted after DIPCH 
in order to form a list of the joints and their coordinates 
for later use (Thsre is no hele), 

6 = LANGIAGE + Icetran 

7 «» PLOW CHART + As shown in Fig. (45.5) Appendix A. 

8 » USAGE + CALL GRID with the following input variables in CGMMON: 


NJ 
MUAOTL 
with the sane neaning 
rUHOIZ 
as in DIPONE. 
HOLES 
VELEN 


9 = COTPUT - Output in COMMCH includes: 
HiCCOR(I) = x-coordinate of joint I 
VECOR(I) = yecoordinate of joint I 
NUJO(I} = Mumber of joint I =f 
1Q— LENGTA - 
Seurce Program « 43 Ieetran instructions 
Object Frogran « 30 Binary cards 
lie CALLS « Hone 
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ige Arr; ICATIONS - SUBSYSTa" PLATES 

13— LONG DESCRIPTION « ALL that it has to be said here is that 
the numbering of the joints rroceeds from the bottom to the 
top and from the left to the right. Although the flow diagram 
is the same es in MUGRID the vrocedure followed differs appre- 
elably. 








1 + NAME - Subroutine GRID 
& @ AUTHOR - John Kolliniatis 
3 = DATE - April 20, 1966 
4% = TYPE - Specific 
§ @ SHORT DESCRIPTION - This subroutine is used in the case of a 
rectangular plate with a hole in order to mmber the joints of 
the equivalent grid and fine their horizontal and vertical 
coordinates, 
6 ~ LANGUAGE ~ Icetran 
7? « FLOWN CHART - As shown in Fig. (495), Appendix A. 
8 - USAGE ~ CALL MUGRID. with the necessary input values in COMMON 
Input variables are: 
HOLENT 
VELENI 
LOLENZ 
VELENZ having the same 
NUHOJ] moaning as in 
RUVESL Subroutine DIPOHO. 
NUEOIZ 
WUVEJ2 
WHiG3 
NUOVEJ3 


{OLEN3 = HOLEN « HOLENI = HCLEN2 
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VELEN3 = VELEN « VELENI - VELEN2 
9 «~ UTHIT - Output in COMMON includes four dynamic arrays of size W 

Where iJ being in COMMON itself, has the same meaning as in 

DIrQdO, These are: 

NUJO(I)= Number of Joint I. 

HOCOR(I) = Horizontal coordinate of joint I. 

VECOR(I) = Vertical coordinate of joint I, 

JORENU(I) + Release nmamber of joint I. 
l0— LENGTH 

Source Program -107 ICETRAN instructions 

Object Program - 45 Machine instructions 
ll» CALLS ~- lione 
l2— APPLICATIONS ~ SUBSISTEM PLATES 
L3— LONG DsSCRIPTION 

This subroutine which is called by the Main program after 
subroutine DIPrOHKO has been executed is used in order to find all the 
nodes of the equivalent grid and their coordinates relative to a 
system of coordinates located at the left lower corner of the plate. 
The mmbering of the joints is as shown in Fig. (C1 ) ives it 
proceeds from the bottom to the tep and from the left to the right, 
in addition if fills with zero’s the elements of the dyna- 

mie array JORENU which is supposed to be filled iater by another 
subreutine by characteristic numbers indicating the kind of the joint. 
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5.6. SUBROUTINE MEIANH 
3 «= NAME «+ Subroutine NEINNE 


2 ~ AUTHOR = John Kolliniatis 

3 + DATE - May 7, 1966 

4% » TYPE = Specific 

5 + SUORT DESCRIPTION + This subroutine is called after DIPONH 
and GRID have been called in order to supply to other subrou- 
tines called later inforsation about the number of any nexber 
amd its mirus and plus joint 

6 « LANGUAGE - Icetran 

7-HLOW CHART - As shown in Fig (45.7) Appendix A, 

8 « USAGE ~ CALL MELNNH with the following variables and arrays 
(defined in previous subroutines). 


HGCOR iJ 
VECOR NUSOMB 
ROJO NUV EMB 
NUME NUREDI 
NULEDI 


9 ~ @TPUT - In COMMCN ineludes the following: 
HUMEMB(I) = Number of member I = I 
MUMLJOCI) = Number of minus joint of member I 
NUPLWJO(I) = * « plus eo a if 
KINDJO(I) = kind of menber I 
KINDJO(I) = 0 side bean 
KINDJO{I) = 1 diagonal 





10° LE.UTH - 
source Program - 71 Icetran tructions 
Object Program - 53 Binary cards 

lle CALLS ~ Lone 

A2— APPLICATIONS — SUBSYSTEM PLATES 

13— LONG DYSCRIPTION + The subroutine exacines whether or not a 
joint is eligible for being a mims joint of a member, If it 
ds the mmber of the corresponding plus-joint 3g found and 
stored together with information indicating the kind of the 
menber. 
Although the flow diagram for this subroutine is the same 
(macroscopically) with the one for MEINHO, the detailed rro- 
cedure is different. That was the reason for writing two dif- 


ferent subroutines. 
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5.7. SUBROUTINE SRINHO 
1 - NAME ~ Subroutine SELNEO 
@@ AUTHOR - John Kolliniatis 
3 = DATE - May 1, 1966 
4 » TYPE - Specific 
5 @ Ssort DESCRIPTION + This subroutine is used in the case of a 
rectangular plate with a hole. Its task is to assign sn identi- 
fication mimber on each member. Find its plus and sims joint 
and store a characteristic number (side beam or diavonal). 
This information is essential for generating the necessary 
information for STRUDL (MEMBER INCIDENCES) as woll as for 
Supplying data to subroutine MEWPRO for the curvose of compe 
ting the properties of the members, 
6 + LANGUAGE «+ Icetran 
7 - FLOW CHART ~- As shown in Fig. (A5.7), Appendix A. 
§ =~ USAGE u CALL MEINHO, with the necessary input values in 
COMMON, Input variables are: 
RUJ OJ1 
1052 
WAQ3 
NOViJL having the sane 
NUVEI2 meaning as in 
NUVES3 SUBROUTINE DIPOHKO, 
RUROMB 
NU VEMB 
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iURIDI 


NUME 

NJ 

HOCOR having the same 
VECOR meaning as in 

NUS O SUBROUTINE MUGRID. 


9 = GTrUT ~ Output in COMMON includes four dynamic arrays of size, 
NUMB, namely: 
NUMEMB(I) = Number of member I = I. 
NUMTJO(I) = Kumber of mims joint cf member I. 
NUMLJO(I) = Number of plus joint of member I. 
KINDME(I) = Kind of member 1. 
ALNDME(I) = 0 Af member I is a sideboam 
KINDME(I) = 1 Af member I is a diagonal 


1O- LENGTH 
Source Program « 139 iCETRAN instructions 
Object Program 85 Machine instructions 


lle CALLS + None 
l2~ APPLICATIONS - SUBSYSTEM PLATES 
13- LONG DESCRIPTION 

The task of this subroutine has been explained in the short 
description above. For the purposes of explaining how this proscram 
functions, the members have been characterized as belonging to 
one of four teams depending upon the direction of their minus-joint 














— SS —_ 
7:2: —— ‘“@—n 
_—— ae : —. 
— - —————— ce el at A © 

a 

ee — me? i 
Ml, a MF * Mh & © 
\ ot -§ em Qe oe eee @ 
= = (ee 

—= lee) (ee 
i — += — Sa 


——aee a * - <li 
ee eel * 
as a 
— “a .- my 


ame = 
wK—l= -.: — —SS = eee ae) eee 


et ee © tees ce oe) oe 
A Ai: § iT eli Aili. 


——=—=—=—=—_ a FS ee Fm ee oe oe 





he 


a* 


te plusejoint vector. These four teams are: 
(4) vertical members. 


(44) horizontal members. 

(444) right diagonal members. 

(4v) left diagonal members. 

The mmbering of the members proceeds from the bottem to 
the top and from left to right in each one of the above teams. 

The sequence of the different teams is as listed above. In other 
words if there are 20 horizontal members and 30 vertical ones, the 
bottom rightenost, right oriented diagonal member will have the 
number 51 and so forth. 

The method that this subroutine follows is to consider 
all joints that are legitimate mims joints for a member in se- 
quenee (joints on the left edge of a hole are not legitimate minus 
joints for tight-diagonal menders except the upper most one for 
examble) and compute the corresponding plus-joint. The mmber of 
the joint as well as the rims and plus joint mambers and the kind 
of the menber information are stored in the apprepriate array 


elements and the program proceeds with the next member. 











lL = NAME - Subroutine FOBOC 

2 - AUTHOR - John Kolliniatis 

3 = DATE - May 2, 1966 

i TYPE - Specifie 

§ = SHORT DESCRIPTION ~ This subroutine is used for the purpose 
of matching the actual supports, with the nearest node of 
the grid. 

6 «= LANGUACE ~ Icetran 

7 = FLO? CHART ~ As shown in Fig. (45.8) Appendix A 

8 ~ USAGE ~ CALL FOEOC, with the necessary input values in COMMON, 


Input variables are: 


HOCOR 

VECOR seme meaning as in previous 
HOO SUBROUTINES . 

MS 


EX = Distance along the x-axis between the center of coordi- 
nates and the particqlar point support. 

WHY = Distance along the Y-axis between the center of coordi- 
nates and the particular point support. 

KIND1= Kind of point support: 
1 - if support is hinged 
2- if support is fixed 

9 = OUTPUT ~ Output in COMMON ineludes 
KINIVO: An array of size WJ containing numbers with the same 
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meaning as above for each node, except that zero inJicates a 
regular node (no releases). This information is necessary for 
subroutine JOINRE in order to produce the proper "JODNIT CCORDI- 
NATES” and "JCLUT RELEASES" statements. 


1Q-] LENGTI 


Source program « 27 ICETRAN instructions 
Object prograa ~ 30 Machine cructions 


ll. CALLS - None 
12 -APPLICATICN ~ SUESYSTEM PLATES 
13~ LONG DESCRIPTION, 

For a given point support the distance from oach joint is 
computed and the nearest joint is located. If this joint is not 
already used as a support, it is considered as such. 

if the joint already is a simple support and the point sup- 
port is fixed, the joint is reclassified as a fixed one. 

After the appropriate joint has been located, information 
is stored in <INDJO regarding its kind, 








§.9. SU as COTIN: 





1 ~ NAME ~ Subroutine LIBOd 

20 AUTHOH - Jonn Kolliniatis 

3 = DATS = May 2, 1966 

4 =» TYPE - Specific 

§ © SiORT D2SCRIPTION - Tais subroutine takes a line boundary 
comlition and identifies it with a set of grid nedes. Only 
lines along the 4 edges of the plate, are considered. Any 
other line boundary condition can be approximated by a set of 
point supports, but the User has to do it. 

6 « LARGUAGE - Ieetran 

@ « FLOX CaAdaT - As shown in Fig. (A5.9) Appendix A. 

8 ~ USAGE ~ CALL LIB0C with the necessary input in COMMON (or 
BQUIVALENCED to other variables in COMMON). Input variables 
are: 

KinOJO 

KOCOR 

VECOR same meaning as in ;O50c. 

NUJO 

iJ 

EXL = Starting (left side) x-coordinate of line boundary 
condition . 

WHY] = Starting yecoordinate of line boundary condition, 

EX2 = Ending (right side) x-coordinate of line boundary 
condition (EQUIVALENCED to EX). 

WHY2 = Ending y~ocoordinete of line boundary condition. 








(BQUIVALENCED to WHY) . 
KIND@ = Kind of line boundary condition (BQUIVALENCED to KLNED) 
KIND2 has the same meaning as KINDL in PORCC, 
9 =~ GUTIUT ~ Output in COMMON includes KINDJO with same meaning 


as in POROC . 

10< LENGTH 
Source Program « 61 Ieetran instructions 
Object Program « 32 Machine instructions 


lie CALLS ~ None 

120 APFLICATIONS ~ SUBSYSTEM PLATES 

43~e LONG UDESCRIPTICN —- LIBOC is to be executed after POBOC if 
there are any line boundary conditions. It examines what joints 
satisfy the ecuation of the line boundary condition and stores 
the value of KIND2 into the corresponding elements of XINDJO(I). 
Of course before it does that it examines whether the line 
boundary condition coincides in all cr in part with one of 
the 4 edces of the plate, If at any tine information is to be 
stored into a particular elenent of KIHINO which already con- 
tains a mmber other than zero the larger value if finally 
kept into the element. In other words if we have two line 
boundary conditions on the sane or different edges touching 
each other at only ene point and one of them is hinged and 
the other is fixed, the common joint will be considefed as 
fixed. 








§o10. SUSROUT F CIMLO 


1 = NAMS — Subroutine FCINLO 





2» AUTHOR - John Kolliniatis 

3 - DATE = May 3, 1966 

& e TIPE - Specific 

§ = SuOkT DESCRIPTION - Initially this subrovtine was written in 
order to take care of the point loads, Later the subroutine 
has been augmented in order to take also care of the displace- 
ments at specified points. 

6 » LANGUAGE - Icetran 

7 ~ FLOW CHART - As shown in Fig. (45.10) Appendix A. 

8 » USAGE - CALL POBCC, with the necessary input variables in 
COMMON. Input variables are: 
NUJO 
NJ With their usual 
LOCOR meaning . 
VECOR 
KIRDJO 
MUS OL 
POLOM1 = Point lead magnitude. 
EXLCAD = xecoord of point of application (SQUIVALENCED TO EX) 
WYLOAD = Yecoord of point of application (BQUIVALENCED TO WHY) 
WOPOLO = Nuaber of point loads 
NLOTY = Type of loading 

i= Z FORCE 








2 = X COMENT 
3 = Y MOMENT 
KIQUT = A variable that is set by CUT and explained there. 


9 - CUTRUT - Output in a punched form consists of a set of “JUINT 
LOADS* “PRINT DISPLACEMENTS" and *PRINT LOADS" statements. 

10. LENGTH 
Seuree Program « 60 ICETRAN instructions 
Object Program - 46 Machine instractions 

die CALLS « NOSUIO 

d2— APILICATIONS — SUESYSTSM PLATES 

13» LONG OSSCRIPTION - The nearest non-support joint is found 
in the same way as in rCB0C. Then according to the value of 
XLOUT the appropriate statements are punched out. The subroue 
tine itself can handle more than one point loads but in 
SUBSYSTEM plates it is used for one lead each tine. 











SUBROUTINE NOSUJO 





 - NAME ~ Subroutine NOSUJO 
2 » AUTHOR »- John Kolliniatis 


~ TIPE « Specific 


un om 
t 


SEORT DESCRIPTION - This subroutine is called after LIECC and 
POBCC have been called in order te generate a set of non-support 
joints 

6 - LANGUAGE - Ieetran 

7 - FLOW CHART - Ag shown in Fig. (A5.11) Appendix A. 

8 = USAGE = CALL NOSUJO with the necessary input variables in 
COMMON Input variables are: 


JO 

HOCOR with the same meaning as in 
VECOR every other subroutine . 
XINDJO 

NJ 


9 = CUTPUT —- Qutput in COMMON includes : 


NGJO1 with the same meaning as 

HOCORL NUJO, HOCOR, YECOR, NJ 

VECORi respectively but for 

NOIZE nonesupport joints only, 
10~ LENGTH 


Source Program « 
Object Program . 
11. CALLS ~ None 


21 Icetran instructions 


25 Machine instructions 








‘A3- The value of KININO is examined for each joint. If it is zero 
| an Sas bosuneasGaneuet yum | 


- - 








5.12. SUBRCUTINE JOIN 





1 ~- NAME - Subroutine JOINRE 

2 « AUTOR - John Kolliniatis 

3 = DATS - May 2, 1966 

4 = TIPE - Specific 

5 = SiiORT DESCRIPTION - Subroutine JOINRE 4s used in order to 
write all the introductory STRUDL statenents, the JOINT 
CCORDINATES table, the JCIUT RELEASES statement and the 
MEMBER INCIDENCES table, 

6 «= LANGUAGE - Icoetran 

7 ~ FLOW CHART - As shown in Fig (A5.12) Appendix A. 

8 «= USAGE ~ CALL JOINRE, with the necessary input variables in 


COMMON, Input variables are: 


WWJO 

HOCOR with the same meaning 
VECOR as in all previous 

NJ SUBROUTINES . 

KINDO 


JURENMS = Joint Release mmber 
9 = CUTPUT ~ Qutput in COMMON includes all the statements montio~ 
ned in 5 above, in a punched form. 


108 LENGTH 
Source Program « 89 ICETRAN instructions 
Object Program - 68 Binary cards 


Lis CALLS - None 
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120 APPLICATIONS = SUBSYSTEM PLATES 
43~ LONG UESCERIPTLION:- The intreductory statanents are punched 
first. Then according to the value of KINDC a table of JOINT 
COORDINATES is punched with the proper EXISTENCE DATA, supe 
port,if joint is a support. Subsequently according to the 
value of JOREM (J) set by LIEOC the proper JOINT RELEASES 
statements are punched. 
In particular if 
JORENU(J) = 0 3 Joint J, releases Moments X,Y 
JORENU(J) = 1 3 Joint J releases Moment Y 
JORENU(J) = 2 3 Joint J releases Moment X 
JOREM (J) = 3; Joint J releases Moment X,Y 
Finally the "MEMBER INCIDENCES" table is punched out. 
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5.13. SUBROUTINE MEXPRO 
1 ~- HAHE = Subroutine HEXPRO 
2 AUTHOR ~ John Xolliniatis 
3 + DATE -~ May 2, 1966 
4% .~ TYPE ~ Specific 
5 @ SORT DESCRIPTION ~- This subroutine computes and punches STRUDL 
statements for the properties of the members of the grid. 
6 » LANGUAGE ~ Icetran 
7 » FLOW CHART ~ As shown in Fig. (A5.13) Appendix A, 
8 » USAG@ - CALL MEMPRO, with the necessary input variables in 
COMMCN, Input variables are: 
CMB 
LUVEMB 
NORTOL Have same meaning as in 
NUME ary other subroutine, 
NUMIJO 
TUPLJO 
NUMEMS 
HOCOR 
VECOR 
PRATIO = Poisson's ratio 
TiICK = Plating thickness 
9 - OUTPUT - Tha output conists of a set of acceptable STRUDL sta~ 


tements that deseribe the properties of the various menbers. 








1O—- LENGTH 
Source Frogram - §3 ICETRAN instructions 
Object Frogran < Je Binary cards 

lie CALLS - None 

12— ArPLICATIONS «- SUBSYSTEM PLATES 

13~ LONG DESCRIPTION + The kind of each member is exariined (side 
beam or diagonal) and then efter computing the length of it 
and all the necessary dimensionless ratio's, the properties 
of the member are computed by using equations 2,18 through 
2.23 of chapter Id. Finally all the above information is 
punched out, forming part of the STRUDL input. 
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5.14, SUBRCUTINE LINELO 

1 =» NAMB + Subroutine LINELO 

2 « AUTUOR - John Kolliniatis 

3 ~ DATE - May 3, 1966 

& - TYPE ~ Specific 

5 © Si:ORT DESCRIPTION » Subroutine LINZ is used in order to find 
the appropriate point (nodes or points along the side members), 
Where concentrated loads have to be applied in order to appro- 
ximate a specified line load. It also handles cases of cutput 
along a straight line. 

6 + LANGUAGE - Icetran 

7 = FLOW CHART - As shown in Pig. (A5.14%) Appendix A. 

8 » USAGE CALL - LINELO with the necessary input variables in 
COMMON. Input variables are: 
HOCCR, VECOR, KINDJO, NUMEMB, NOMEJO, NUPLJO, NUVEMB, KUHOMB, 
KJ with the save meaning as in previous subroutines. In addition: 
XA = Start (left) x-coordinate (BQUIVALENCED to EX) 
YA=s * # —_—-yeecoordinate (EQUIVALENCED to WiiY) 

end (right) x-coordinate (EQUIVALENCED to EX2) 


u 


Xb 
Ye« * 0 y-coordinate (BQUIVALENCED to WHY¥2) 
STALD = Start loading Magnitude (EQUIVALENCED to FOLOM1) 
ENDLO = End loading Magnitude 
NLCTY = Same meaning as in PCINLO 
LTYFE = 6 if uniform load 

1 if linear load 


10- 


ii. 
l2= 


OT UT -« Output in a punched form includes STRUDL statements 
that will be explained in the chapter for LINE since this 


subroutine serforms all I/C operations for LINELO 


LEGGTH 

sources Program « 230 Iecetran instractions 
Object Procram - 284 Binary cards 

CALLS - LINS 


APPLICATIONS + SUBSYSTEM PLATES 

LOKG DESCRIPTION 

If the line load does not coincide with a set of members, then 
the intersections of the line load projection with all menbers 
intersected by it are found. The number of the intersected 
member and the local and global coordinates of it are stored 
for later use by LINE. Special consideration is also given to 
the case where the line passes thru a joint or a small area 
arround it in order to be considered as a JOINT LOAD instead 
of a LINE LOAD. 

In case that the line coincides with a set of members then 
the mumber of the coinciding (with the line) members as well 
as local and global coordinates of the start and end loading 
point for each menber, are stored for later use by LINE. 

In case that this subroutine is used for a line output certain 
parte of it are by-passed. 

Finally subroutine LINE is called which in addition to other 


tasks performs also all punching necessary for LINELO, 
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5,15. SUBROUTINE LINE 


4 wm HAME ~ Subroutine LINE 





@= AUTHOR ~ John Lolliniatis 

3 ~ DATE - May 3, 1966 

4& - TYPE ~ Specific 

5 = SUORT DESCRIPTION - Subroutine LINES is a logical extension of 
LINELO separatedfrom it in order to get arround the difficulty 
of the precempller to accept a larze mumber of symbols. It is 
ealled by LINELC in order to coapite the different grid loads 
that are ajyrivalent to the plate line loads specified. In addi- 
tion it performs all the 1/0 functions for LINSLO for both 
the case cf line loads ard line outputs. 

6 » LANGUAGE - Icetran 

7? - FLOW CHART - As showmm in Fig. (A5.15) Appendix A. 

8 ~ USAGE ~ CALL LINE (KINDLO, AGLLD, YGLLD, NULOME, XLCCAL, 
YLCCAL) 
in addition all tre input in COMMON tn LINELO must also be in 
COMMON in LING plus tho variable KIOUT. 
The values transmitted between LINELO and LINE as argunents 
have different meanings according to the value of the variable 
KINDLO, £INOLO itself is an array of size NULDS which contains 
munbers 0,1,2, selected as it was indicated in the flow dia- 
gram for LINELG. (See Appendix EB for meaning of NULDS) 
In particular: 
(4) If KINDLO(J) = 0 (Member concentrated load) 





3 


-” 


XGLLD = Glicbal x coordinate of concentrated member load 
YGLLD = Glehal y coordinate cf : ° » 
NULOVE = Mamber of loaded member 
ALCCAL = Local x coordinate of concentrated member load 
YLOCAL = Local y 2 * ¥ ® ” 
(44) If KINDLO(J) = 1 (joint load) 
XSLED 
= Global coordinates of loaded joint 
YGLLD 
NULOME = Number of loaded joint 
(134) If <INDLO(J) = 2 (Member uniform or linear load) 
XGLLD = Clobal start of mexber line load 
YGlID= ona * " C 
XLOCAL = Leeal start ef member line load 
YLOCAL = Local end * » », & 
HULOME = Number of loaded member 
in this case it is necessary to say that XGLLD, XYGLLD, 
XLOCAL, YLCCAL, all refer to the changing ccordinate, 
singe in this case either the x or y coordinate remains 
unekanged. 

9 = CUTPUT ~ Cutput in a punched form depends upon the values of 
different key parameters, Table (515.1) gives the different 
combinations of these parameters 

10—- LENGTS - 

Source Program - — 269 Icetran instructions 


Object Program «- 250 Binary cards. 








APPLICATIONS — SVESUSTEN PLATES 

LONG DESCRIPTION -~ If KINDLO(J) = 2 the list of the points that 
will be used to approximate the loading is shorted according 
to the distance from the left end (or lower end if load line 
is parallel to y axis), 

If KINDLC(J) = 2 this step is unnecessary end it is skipped. 
Zach ene of the above points takes so mich of the Joad as it 
is applied batween the middle points ef the segments joining 
it with its neighbouring points. Of course this is not the 
ease when KISDLO(é) = 2, for all J's, In this case each member 
is loaded with the exact linear or uniform load acting over 
it, Subsequently the necessary output is punched acecrding 


to the values of the contrel variables as in table (5.15.1). 
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5.16, SULBCOTINE ARDALO 
ecg 


i. 
2-= 
3- 
4 
5 = 


6 - 
7 - 
8 - 








Ali = ouoreuline AnLALO 

AJDT = Jolin Agliiniatis 

UATE = hay 3, 1966 

Titi © specifie 

wulal vesCRirTlON - Subroutine AniaLo is used in order to treat 

area loads. tne metned is used to divide tne Loading among 

the members covered by the area load and then call SUDAGUTINE 

LIABLU. Unly loads with perimeter rectangular with the edges 

parallel to tne edges of the plate can be treated, 

LANGUALUL - Icetran 

FLGW ChaRT ~ As snown in Pig. (45.16) Appendix A. 

USAGE - CALL ANZALO with the necessary input variables in 

CUiNON. Input variables are: 

HOCGR, VECOR, NUMIWJO, KUFLJO, MUVENS, SUNOHB, NLOTY, with 

their usual meaning. 

4A = x-coordinate of lower left corner of area load 
(EQUIVALENCED to EX) 

XA = yeocordinate of lower left corner of area load 
(EQUILVALENCED to WiY) 

A5 = xecoordinate of upper right corner of area load 
(EQUIVALENCED to EX2) 

YR = yecoordinate of upper right corner of area load 


(EQUIVALENCED to WhY¥2) 
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LOMSTA = Starting loading magnitude (at lower or left odze) 
LTZPE = A variable that indicates the kind of the load;in 


particular: 
LTYPE = 0 Unifor: load 
LTYPE = 1 Load linear along the xeaxis 
LTYFE = 2 : g " y-axis 
10- LENGTH 
Source Program « 130 Icetran instructions 
Object Program - 109 binary cards 


11 CALLS - Subroutine LINELO 

12— APFLICATIONS — SUBSYSTEM PLATES 

13e LONG DESCRIFTION «+ The members covered by the load are first 
located. If LIvPi is 0 or 1 the load is distributed axong 
horizontal members, if LTYFE = 2 the load is distributed axnong 
vertical members. Then subroutine LINELO is called in order 


to punch the necessary output statenents. 
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5017. SUBROUTINE CUT 


1 « 


NAME «- Subroutine CUT 


2 ~ AUTHOR » John Kolliniatis 


3- 


h « 


5 « 


DATE - May 3, 1966 

TYPE - Specific 

SHORT DESCRIFTION - This subroutine examines the value of 
KIOUT, which is set by the usor’s data or COMMAND (4f CDL 

4s used) and according to its value the necessary steps are 
performed by GUT or other subroutines in erder to produce the 
specified OUTPUT. 


6 ~ LANGUAGE <- Ieetran 


ia 


FLOW CHART « As shown in Fig. (A5.17) Appendix A. 


8 » USAGE ~ Cal] FOBOC with the necessary input variables in 


COMMON. Input variables are: 

KiCUT imdicates desired output. In particular ifs 
KIOUT = 2 - LIST REACTIONS 

KIOUT = 3 ~ LIST POINT LOADS XA.. YA.. 

KIOUT = 4& = LIST LINE LOADS XA... YAese XB... YB 
EIGUT = 5 ~ LIST LOADS 

KIOUT = 6 — LIST PCINT DISPLACEMENTS XA... YA 
KICYT = 7 = LIST LINE DISPLACEMENTS XA... YAeos XB... YB 
KICGUT = & = LIST DISFLACSMENTS 

OUTZA = x ard y coordinates of point 
CUTYA = where output is requested 
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Bo 


coTxa | 

| x and y, start and erd coordinates 
CUTYA | 

\ of line along which outocut is desired 
CUTXB 

| 
QTY | 


ae 


9 ~ WTFUT - The output is in a form of punched STRULL statements 
that take care of any desired collective or selective output 
that the USER desires. 


1Q— LENGTH 
Source rrogram - uO Icatran instractions 
Object program -« 27 Binary cards 


lie CALLS - LINSLO, POINLO 
42— APPLICATIONS - SUBSYSTEM PLATES 
13@ LONG DESCRIPTION - Flow diagram and short description indicate 


completely now this subroutine functions. 
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CHAPTER VI 





Several problems have been runned, for the purpose of deb- 
buging individual subroutines and examining the interaction of then. 
The resuits have been excellent regarding the cenerated output. 
Gniy the input and output for a very simple croblem wes included 
in Apperxliz *, It covers the case of a square plate with two clamped 
and two hinged edges under linear, lateral, area load. it has teen 
selected since is one of the solved problems in (11). 

There is nothing that can be said on the overall accuracy 
of the method since the nonwavallability of STRUDL prevented the 
author from solving the problem completely. However judging from 
the results cf (1), where a grid has been generated mamally, the 
larger deviation frem values measured by experiment are of the 
order of 15%. 

After the disclacements have been found the moments and 
shears can also be evaluated (more correctly) from the following 
relations given in (2). 


Mx =0 ( Sonate) (6.1) 
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= -0( +b 3 ) (6,2) 
namie: os ‘ hop) (643) 
SR oa oe ) (6.4) 
whore 
w = lateral displacenent 
Diceeele Le (6.5) 





1QCi- BD) 


Before we state some recommendations for future expansion 
of this work, there are several other points that mst be said 
regarding its completion, First the output from this program has 
to be tested in STRUDL when it becomes available. Second the CEL 
program (a first draft of which is shown in Appendix D) has to be 
compiled, tested with dummy subroutines and the subroutines of this 
program. This will make the srogram for rectangular plates in fle- 
mre an integrated one. 


In the very restricted category of rectangular plates the 
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bending ercblem is not the cnly problem, There also axial loads, 
instability cr combinations of them. A further expansion of this 


work, therefore, would be to write sinilar programs for these cases 
and then put everything together, It would really be a great contri- 
bution to engineering to be able to analyze any plate with any 
combination of lateral and axial leads in a matter of minutes. 
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CALL 
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=z 
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CALL 
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MAIN PROGRAM 
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POINLO 





756 















ARE THER 
kN Y 


LINE LOADL, 










ARE THERE 


ANY 
AREaArconns * 


READ 
REQVESTED 
OvTPuUT 


ENDO 


Bien aoe 


READ wtine 
LOAD 
DATA 


READ AREA 
LOAD 
DATA 


CALL 
ARE ALO 





cont. 


(Ae 





INCREASE 
NuHoyr? OANUVeY) 
al 


Yes 








DECLARE DYN. 
ARRAY 


COMYUTE FRARCTIONA 
Ne 
OF DIVICion POINTS 
TRUNCATE 











POssiBLE to 
INCREASE EITHER 
'NuHost' or 

"Nuovete 





COMPUT 24 Uae 
oO 
DYN. ARRAYS 






DEFINE GLOBAL 
DYN. ARRAYS 





Fige (A5.2) 


SUBROUTINE DIPONH 









Vile 










COMPUTE 
NEW Manuyso’ 


DECLARE DYN. 
ARRAYS 


COMPUTE FRACTIONAL 
No 
OF Division POINTS 
TRUNCATE 
COMPUTE SIZES 
OF GLAL DywNw ARRAYS 


PEFINE ALL GLOBAL 
DYN. ARRAYS 


RETURN 


Fig. (A5. 


SUBROUTINE _DIPOHO 


(ier 





HOcoR CI)= 


- WOCOR (1-4) 





Iz ivi 


ENTRY 


DECLARE 
DYN. ARRAYS 


COMPUTE LENGTH 
OF 


SIDE RPEAMS 


LOCATE S$uB-PANEL 
WHERE Iyer is 


SET NEW HOcoAR(:2) 


SET MEW VECOR CE) 


SET 
Nuzoczr)-T 


i$ JOIN 
NuUTO (C2) On 
THE UPPER 
Epceé? 


oc S 















IS 
JoinT T 


ON THe VePER 
RIGnhT 


Conners 


YES 


SET FOR ENVCL) BRU 
TO ZERO FORALL I'S 


SUBROUTINES GRID AND MUGRID 


VE COR (r41)=2 0 


T= 1* 1 


79. 





80. 


ENTRY 
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DECLARE OYNAMIC 
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DECLARE DYN, 
ARRAYS 


DEFINE LOCAL 
DYN. ARRAY NV 





INTRO OVCTOR 
‘$TANOL' 
STATEMENTS 











4-0 aAWwP 
| STaatwity s**t rar, 









Join 
List 
Ext Austen? 


Picwu-ve 
NEXT FOINT 








YES 


7 4 


WRI TE 
APPROPRIATE 
‘JOINT RELEASES 








Fi A5,12 


SUBROUTINE JOINRE 





START with It mage 





€$ 
ExUAUSTED 





WRITE 


NUMEMB 


NYUM\Jo 
Nv@L_LJo 








DESTROY NYU 


Wal A5.12) cont. 


SUBROUTINE JOLNRE 


87. 





DECLARE OV N. 
ARR ANS 


DEFINE LOCAL 
DYN ARRAYS 


ZEAGQ IN ALL ELEMENTS 
OF REINER & TOINER 


CONSIDER FIRGT 
GAID UNIT 


GOMPYTE REQUIRED 
PROPERTIES 


Apd fRotERTIES OF 
Six MBAS OF UNIT TO 
CONTENTS of REINER, TOINER 


E> 


Yes 










WRITE 
NECESSARY 
STATEMENTS 





DE *TROY ALL 
LOCAL DYN. ARRAYS 


Fic. (A5.1 


SUBROUTINE MEMPRO 


CONSIDER 
NEXT 





ON AT 


88. 





WES 


mo 
COincidve witn 


AWY MBAS 
























FIND AND 
STORES 
LOCAL AND 


GLo8 Av 
CooRnDInATes 


Finp 
IN TEQSE CT's 
with VERT. 
MARS 
STOAE LOCAL 
AND GLOBAL 
COOK O:wates 


OF LoApeo ABR 


ENTRY 


DECLARE DYN, 
ARRAYS 


DEFINE 
DYN. 


LOCAL 
ARRAY 5 














1s 
LINE 
PARALLEL | 
TO As EDGE: 


NO 


ON THE LINE 


STORE NHR 
OF so;init 
41 NoLo(s)= Ih 


Yes 













Fino AND 
STORE 
LOCAL 

AWD GLOBAL 
COORDINATES 
oF LOADED K8C 













INTERSECTIOMS 
WITH HORI 2, 
MBS 

STORE LOCAL 
AND GLOBAL 
COO ADIMNATES 













MINPLO(I)= 2 


RECEAT Untr 


Avi POINTS 
FOUND 


DESTROY 
LocAt ARRAY 
RETLRAN 


Fig, (A5.14) 


SUBRCUTINE LINELO 
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Fino 
LNOECASECT’Y 
Witia VERY 
MEMAGR)D 
ST#AE BC 
Ane GroBdAr 
COOK ONAF. 








Find 
IMTEASEcTNS, 
with WOR. 
MEMBERS 
$TORE tCLe 
AAO GLOBAL 


COCRDiM:»4 


141 NOLO (13) 
ao] 











DECLARE DBYN 
ARRAY 5 


SHORT LisT 
oF POINTS 


YES No 


COMPUTE 
LOADINGS 


EXAMINE VALVES 
OF 







Win DLO 
NLOTY 
LTYE 
KiouT 













APPAOPRIATE 
STATEMENTS 


AS IN TABLE (4.1) 


SUBROUTINE LINE 








PDISTAIBYVTE LOAD 
AMONG 





FinDdD; wWORIE. 
COVERED BY LOAD 


WORE MARS 


MaRS 


ENTRY 


DECLARE DYN 
ARRAYS 


$ 


DiSTRIQvTE wroavrpn 
AMONG VERT. MBRS 










FIND VERT, MERS 
GOVERED BY LOAD 


SET NECESSARY 
VARIABLES 


= 


Big. (AS. 


SUBROUTINE AREALO 


91. 





92. 


DECLARE OYA. 





E $ 


WRITE 
Yt RODucTro Ry 
STATEMENTS 





"(QINT 
REActTions’ 





PRINT 
DisPLACE Meme? 
LOADS 


Fi A5.1 


SUBROUTINE OUT 





= 
JJOB 7 KOMI NAT 1S —— 





/RESTORE OLD 69 
MALI BRARY_ICESI 1B. 
/DUMP 
___ **026 —-_ i 
DYNAMIC ARRAY HOCOR»sVECOR sNUJOsJORENU s K INDJO»s NUPL JO» NUMI JO sNUMEMB 
DYNAMIC ARRAY KINDMEsNUJO1]sHOCOR]>VECOR) SSS 
REAL LDOMSTAsLDMEND 
COMMON HOLENs VELEN» THICK sMANUJOs YOUNG »sPRATIOsHOLEN] sHOLEN2 s VELEN] 
1» VELEN2 9EX sWHY sKIND1 9EX29WHY 29PCLOM] sKILILOsLTYPEsNUHOJ1 sNUVEJ1 9 
2NJ sHOCOR(P) sVECOR(P) sNUJO(P ) » JORENU(P) sHOLEN3 »s VELEN3 sNUHOJ2 s NUHOJ3 
3 »NUVEJ2 sNUVEU3 sNUPLJO(P) sNUM]JO(P) sNUMEMB(P) sKINDME(P) » NUME sNUVEMB 
| 4» NUHOMB »sNURID I »NULEDI sNUPOLOsNUJO1(P) sHOCORI(P) sVECORI(P) »sNSIZE + 
5NULDSsENDLDsLDMSTA sLDMEND sKINDJO(P) sKIOUT sNLOTY 
| FQUIVALENCE CEXSEX1]1sEXLOAD2XAs0UTXA) » (WHY sWHY]9WYLOADS YASOUTYA) 9 | 
ReTND 1s KIND?) s°CEX 29XB » OUTX Ba » (WHY2 > YBsOUTYB) » (STALDsPOLOM1 } 
mee CSCC CF OORMAT (5 F103) 
3 FORMAT (16) 
eee AN: (lates G1 O03) ee ee 
: 5 FORMAT (214) 
6 FORMAT (1494F10.3)-. 
| 7 FORMAT (314) 
. BmORMA Tl 155 3F.10 03g 155 3F1053 ) —— mf 
| 9 FORMAT (1594F 10.3) 
Pp GALL WNIT(QQDUB) 
READ (5592) HOLENs VELENs THICK s YOUNG sPRATIO 
WRITE(692) HOLENs VELENs THICK» YOUNG sPRATIO 
READ (533) MANUJO 
ee 6.» 3) MANUJO. — —— 
READ (594) NHOLsHOLENI» VELEN] sHOLEN2 sVELEN2 
WRITE(694) NHOL»sHOLEN] » VELEN] sHOLEN2 s VELEN? 
. Pay TO (50:60) »NHOL 
























































O CONTINUE _ 
q THERE IS NO HOLE 
HOLEN1=HOLEN 
VELENL=VELEN 
CALL DIPONH 
10 CONTINUE 

















11 CONTINUE 
CALL MEINNH~ 
| 12 CONTINUE 
. fe TO 70 
| 60 CONTINUE 
mere TSA HOE. 














| CALL DIPOHO 
| 13. CONTINUE 
| CALL MUGRID 








ONT INUE 
CALL MEINHO 





' 7O CONTINUE 
: READ (595) NUPSUsNULSU 
WRITE(635) NUPSUsNULSU 
IF (NUPSU) 85985580 
| 80 CONTINUE 
DO 85 J=1+»NUPSU 
READ (596) KIND19EX»WHY 
WRITE (696) KINDI2EXsWHY 
____NUPSUP=1 | 
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56 Nh 


_ CALL POBOG —_ 
85 CONTINUB 
, L:Peeel NI eS) eB ely5.0.9.0. 
50 Peon {NUE 
YON Bere ee OTS, 
READ (536) KIND2sEX1 »WHY1sEX2sWHY2 
ee WALL Feel Coed IND 2a X 1 »s WHY lo BX? > WHY? 
NULSUP=]1 
= CALL LIBOC | 
ioc ONT I NUE 
- SA JOTNRE 
16 CONTINUE 


950 


fee CALL MFMPROO a 


100 CONTINUE 
— READ (557) NUPOLsNULILsNUARL 
WRITE(6s7) NUPOLsNULIL »NUARL 
7 IF (NUPOL) 11051103105 











HeO> CONTINUE 
e110 T=JwNUPOle 
rOouT=1] 
meeUPOLO=)__. oe 7 . tan a 
READ (539) NLOTYsEXLOADsWYLOADsPOLOM] 
WRITE (639) NLOTYsEXLOADsWYLOADsPOLOM1 _ 
CALL POINLO 
mo mGONT TINUE 
Mme (NULIL) 12031203115 
fee CONLINUE 
PeoT20 T=TsNUCTL 
eee KTOUTL=) ————— 
READ (598) NLOTYsXAsYAsSTALDsLIYPEsXBsYBsENDLD 
WRITE(698) NLOTYsXAsYAsSTALDsLTYPEsXBsYBsENDLD 
meee |L_INELO 
meee wen INUE — wH.  _ 
IF (NUARL) 13051309125 
lez SOON TINUE 






































DO 130 T=1sNUARL 
ano Ua. = 

READ (538) NLOTYsXAsYAsLDMSTAsLTYPEsXBsYBsLDMEND 
WRITE(638) NLOTYsXAsYAsLDMSTAsLTYPEsXBsYBsLDMEND 
eee AREALO 
eee ONT 1 NUE 

















READ (535) NUOUT 





WRITE (635) NUOUT 


DO 140 [=] sNUOUT 
READ (539) KIOUTsOUTXAs0OUTYA:s OUTXB »sOUTYB 
WRITE(699) KIOUTsOUTXAsOUTYA sOUTXB sOUTYB 
CALL OUT 




















140 CONTINUE 
141 FORMAT (T85s8H FINISH ) 
WRITE (79141) 
DESTROY HOCOR 
DESTROY VECOR 
DESTROY NUJO __ 
DESTROY NUMIJO 
DESTROY NUPLJO 
DESTROY NUMEMB 
i __DESTROY KINDME 
DESTROY KINDJO 
DESTROY JORENU 
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iro coamaind ‘tel 


ee fe ee ae . 
yee f.>8- -seay 
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Baie : 











ee DES TROYRNU IO}, 
DES Ts OVRICGery 
DESTROY GOR. — = — : —_ 
END 
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e*UZO 
SUBRTGOT INE WD I RUN 
DYNAP IC ARKAY RCCUR sVE CUR ofdifUC »s JORENG see WwOJUC ofthe LIU s wh 1 UO sito 
DYNAMIC ARRAY KINOME sNUJO1 sHUCURI » VECURKL 
COMMON HOLENs VELEN» THICK »sMANUJU 5 YOMLING s PRAT IU sHOLEN 1 s HOLE WN? s VEL EMT 
1»sVELEN2 sEX ®WHY sKINDI] sEX2 9 HY2 sPOLUP 1 sKILILUSLTYPE sWUHUI] swUVEJI.- 
2NJ sHOCOR(P) sVECOR(P) s\UJO(P) sUURENU(R?) sHULEN3 5 VELEN3 sNUHUJZ -NUHOJ > 
3 eNUVEL 2 sNUVEI3 sWUPLIUIP) sNUAITJU(P) sUd4EMB (CP? sKTiUint (RP) so NUME oe NUVLMO 
4 sNUHOMBsNURIDI sNULELI»sNUPULUsNUGJUOL1 (PR! sHUCUKI1 (FH). VECURI(P)NSIZEs 
2NULDS sENDLDsLOviSTAsLDMENDsKINDIU(P?) sK [OUT «NLUTY 
Meet PiKATLO<le/se0 2szat 
1 WRITE (624) 
Seal of XT T 
CONTINJE 
PeeeeRMAT (41H POISSON'S RATIV LARGE sANAEOGY Nerval) 
RATIO=VELENI/HOCLEN] 
ANUJO=MANUJO 
AUVEJI=EASCRT)ANUJO#RATI UX 
AUHOJIEZEAUVEIJL1/RATIO 
NUVEJIL=AUVEJI1 
NUHOJL=AUHOJ1 
mee VER=(NUVEJ1+1)*NDAOJ1 
INCHOR=NUVEJ1* (NUHOJ141) 
meer oMANUJO—-INCVER) lds 5210 
Pee ONT INUE 
<a MAX NUMBER OF JOINTS ACHEIVED 
NUVEJIENUVEJ1+1 
mo TO bu 
mem COUNTINUE ? 
C meer INE POSSIBILITY CF INCREASING HURTZ.IGENTS 
Meet MANUJO-INCHOR) 25e2v53U 
2U CONTINUE 
‘¢ Pee NUMBERK OF JOINTS ACHEIVED 
NUHOJ1L=ENUHOQJ1+1 
BOOTO Sev 
meee ON TINUE 
C MmeerOScIblE TO INCREASE EBIlLIBER [HE HORIZ OR VERT Sos 
mo 10 50 
PemeGONTINUE 
@ EITHER THE HORIZ OR VERT JOINTS CAN BE INCREASEDeMAKE LARGER CHANGE 
Mee (NCHOR-INCVER) 40935235 
35, NUHOJ1=NUHOJ1+1 
Mo, 1O 5u 
40° NUVE J1=NUVEJ1+1 
memeoONT I NUE 
NUHOMB= (NUHOJ1-~1) *NUVEJI1 
Pe veEmMb= (iwUVEJI—-1)*NUHOJI 
NURIDI=(NUHOJ1—-1)*(NJUVEJI-1) 
NULEDI=NURIDI 
NUME =NUHOMB4+NUVEMB4+NUR IDI4FNULLDI 
NJ=NUVEJLENUHOJI 
MePINE KINDME sNUMEsF ULL=HIGH 
DEFINE NUMEMBsNUME »5FULL»sHIGH 
Per ine NUMIJOC-NUMEsFULL> HIGH 
DEFINE NUPLJO>NUME.FULL 2HIGH 


ps! 
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DEFINE 
Viet IN 
EE Nite 
DE EL NE 
betel Nie 
OE te aie 
DEE lite 


wer (NE 
MA 
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TURN 





NUJO= 


VEEURiasNJsF ULL shI GH 
mOCWRI s Nis FULL oe 1GH 
NUJOL “Ss Nols F Uleles ia es 
KINDJOsX\JsFULL>HIGH 
JORENUsNJ»sFULL > HIGH 
NiO sIN s Fi Pigs eile 4 

VECCR »NJsFULL»HIGH 


HOCOR »NJsFULLsHIGH 
” 
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Ge 


SUBROUI TL Wt “lr URS 

DYNAMIC ARRAY FPOCVER s,s VECUR, WOU - OKANO 1 WO INURL 6 RIT JU 5 WU 
DYNA‘ IC ARRAY sTNDME sNUJOI] sHOCUR] = VECURTL 

COMMON HOLENsVELENs THICK sHlANUJu s YUUNG s PRAT IV sHULEWL sHULEWe sVELENI 
LsVELEN2 ste xX sWHY -KINDI ot X2 eWHY 2 -PUOLWa sk 1 Libs Yr meer eee ee 
2NJIsHCCUR(P) »VECUR(P) sJUJU( HK) s JUROR WP eriCLEN Ss VELENSSNURPUUe Meese 
3esNIUVEI2 sNUVES 3 se NUPLIVU(P) eNUMIL UU (HK) ewworehip (PP) sn TINDMIL (RP? so NUME SNUVEMD 
4eNUHOMBs NUIK IDI sWJLEDI] »>NUPULUsNUJUL (FF? sHUCUKI (fF), VecCUKI(P)sNSI Ze; 
SNULDSSENDLD>LDMSTA.LONENDSRINUVIUUCP) sora] ly 

DaMieN ol OND (6 ) BY Go) sion eo 

REAL OD 

lALS.SUBROYTINEA tS CALLED oY FUGK Ite IT GIVES THe NUMBE =o Jo ee 
mor EACH PART OF THE PLATE IF THE AA PCE. Oe voto to 

Pee clr TEDAWIHERE 1S A aOke 

meme RAT TCS w/ Ae) os ae 

WRITE (634) 

meek EXIT 

CONTINUE 

PeRMAT (41H POISSON'S RATIO LARGE SANALCOGCY WOT VAT oe 
HOLENZ=HOLEN-HROLENL-HOLENZ 

PeeeENS=VELEN~VELEN]—VELENZ 

mol) =HOLENIL 

D(2)=HOLEN2 

PS ) =HOLEN3 

D(4)=VELENL 

D(5)=VELEN2 

B(6)=VELEN3 

DO 5 K=1;36 

fee y=D(K)V/SD( 1) 

CONTINUE 
ALPHA=A(4)4+A(5)4+A(6)4A(2)*A(G4IFA02Z) ®A(6)4FA03) FAG) FAC ZI FAL OTS 
me 3)*A(0) 

BET A=letA(3)4+A(4GI4+A(6) 

ANUJO=MANUJC 

me Sc K=-1:6 

Pe =4*ALPHA*((A(1)/A( 8) )*#2) 

IF (ANUJO-ANUJ) 75896 

ANUJO=ANUJ 

MANUJU=ANUJO 

MONTINUE 

mo 10 K=156 

AN) 1*=BETA/ALPHA eSQRT)) ) BETAS ALPHA***¥ 2% eANUJU/SALPHA* 

AN( KI =AN( 1) #A(K) 

N(K)=AN(K) 

GONTINUE 

ee C1) TN G2 IAN 3) 2 

N2Z2=N(4)4N(5)4+N(6)-2 

NIJENLEN2Z-N( 2) 4N (30) 42% (NO 2) 4N(9))-4 

NUHOJ1=N(1? 

NUHOJ2=N(2) 

NUHOJ3=N(3) 

MUVEJLHEN(4) 

NUVEJ2=N(5) 

NUVEJS=N(6) 








————— 
EE _ ~ 

= 
tS cm rm ee 


es 


—_——_ ~~ 














NUHUIQSENUNC OL +iNGHUJ 24+ JH. vse 
NUVEJU=MEIVEDEANUVE JZFNUVEU 3] 

NUVEMbD=(NUVEJURK1) %#NUnUIUR—(WUVLS era) ® GWU ewe) 
NUHOPFib = (NUHOUJO-1)* (wOUVEJOUI-(NUVL JZ See (WUTC le - 7 
WURIDI= (NUVEJU-1)* (NUHOJU=1)— (ROS ee ee ee 


NULEDI=NURIDI 
NUME =NJHOMbt+NUVEMib+NUR TDI +yUbcu! 
ae ens FULLeshIG 


























101. 


a *+xO76 — 
SUBROUTINE GRID 
DYNAMIC ARRAY HOCOR »VECOR sNUJOs JORENU sKINDJO sNUPLJO»s NUMI JO »>NUMEMB — 
DYNAMIC ARRAY KINDMEsNUJO1sHOCOR1 »VECOR] 

COMMON HOLENsVELENs THICK sMANUJOs YOUNG sPRAT LOsHOLEN] »HOLEN2 sVELENI 
L»>VELEN2 sEXsWHY sKIND19EX22WHY2 sPOLOM] sKILILOsLTYPE sNUHOJ] sNUVEJ1 5 
es CN I SHOCOR(P) eVECOR(P) sNUJO(P) « JORENU(P ) sHOLEN3sVEI ENS NIIHO I Qe OD) 3 
3 sNUVEJ2 sNUVEJ3 sNUPLJO(P) »sNUMI JOC P) sNUMEMB(P) sKINDME(P) »NUME »sNUVEMB 
4sNUHOMB»NURIDI sNULEDI »NUPOLO»sNUJO] (P) sHOCORI] (P) sVECORI(P) sNSIZEs . 
SNULDSsENDLDsLOMSTAsLDMENDsKINDJO(P) sKIOUT sNLOTY 
z FOUDVALENCE (NUHOJ bsNUHOUO)-s (NUVE BI sNUVE DO) es a 
a e#® THIS SUBROUTINE GENERATES A GRIDsGIVEN A RECTe PLATE WITHOUT A HOLE 
= CL) Se 
NHOSPA=NUHOJO-] 

— HOSPAN=NHOSPA_ 

VESPAN=NVESPA 

_____sXSPAC=HOLEN/HOSPAN 
YSPAC=VELEN/VESPAN 

J EE 

NVESPA=NVESPA+1 

NHOSPA=NHOSPA+] : 

DO 30 IO0=1sNHOSPA 

I=[0-1] 

DO 20 JO=lsNVESPA 

J=JO- 

R= PANUVEJO+J+1 

Al=I a 

IF (I-NUHOJO41) 12310512 

10 HOCOR(K)=HOLEN a — ee 

Geo 10 13 
MereereNE as 
HOCOR(K)=AI*®XSPAC 
13. CONTINUE 
AJ=J 
Ree CO—NUVEJO+)) 15514215 _ 
14 VECOR(K)=VELEN 

eae GO 10 16 

15 CONTINUE 

pen (kK) =A LWMEYSPAC 

16 CONTINUE 

BOK.) =K 

PoeeONTINUE 

i 
30 CONTINUE 
DO 31 J=19NJ 
JORENU( J) =0 
31 CONTINUE 7 me I — “oS 
RETURN 
END 


me 
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SUB OMT Dif 4b 1 Ned 

DYNAWBIC AKKAY MhUCUR s VECUR sinuUUY sUURLIN aN DIU YS te HOURLY eNO] OU of OIE pe 
DYNAMIC ARRAY <INDwe st§¥UJUL siUCUKI > VECURL 

CUMmmiGyn HYULENs VELEN s THI CReMtANUJU s YUNG s PRAT LUSHULEW 1. > HUE 3 oe 
L»sVELEN2Z sEXeWHYs KIND] SEX! sWHY2 sPULUP Len TLILUSLIYPE snUHUs De N Uy Eee 
2NIsHOCUR(P?) sVECOK(P) sVUJU(P) sJUIENUC Ed sRULEN3 2 YELCIN3 sNJHUIZ UU IS 
3 sNUVES2 sNUVEI3 sNUPLIU(P) swUNLUU(P) siwuOMEMo (RP? eK INOML (CP) eNUrE. NUVEMO 
4 sNJHOMBeNURIDI sitULEDI sNUPULUsNUJUL (FP) sHUCURKI1 (FH) s VECURI](P).NSIZEs 
SNVLDSsENDLD »L DSTA sLOMEND sKINDJO(P) -KIUUT-NLUTY 

E@UIVALENCE (NUHOJ] .AUA@US) .CNUVEDT uve) 

DO 4 J=1.sNUME 

SNOT € | jay 

CON RENUE 

NUJENUHUVIOXNUVEJUH=] 

K=U 

DO 2c JT=leNUJ 

ie wri (HOC PR 141) Osea 

CONTINUE 

J=I-K 

mem) JO( J) =NYvJot]) 

moe lL JOC J) =RUJC(I+1) 

NUMEMB(J)=J 

BO e1TO 2u 

CONTINUE 

K=K+1 

SONTINVE 

HORIZe“EMBERS 

Kee 

NUJ=NUVEJO* (NUHOJU-1 ) 

DO 4C J=1.NUJ 

J=1+(NJVEJO=-1)*NUHOUO 

pert JOC J) =NUJO(]) 

K=I+NUVEJO 

eee ea) ) =NUIO(K) 

NUMEMB(J)=J 

CONTINUE 

DIAGONAL MEMbDERS 

NUJ=NUVEJO* (NUHOJO-1)-1 

K=U 

DO 6v JT=l1 5,NUJ 

me) 1TO=VECOR( I) /VeECORNUveEsD) 

IF (RATIO-“-1le) 455350345 

SON TINUE 

PraeReE IS A DIAGONAL 

J=1+(NUVEJO-1) *NUHOJU+F (NUHUJU-1) FNUVE JUNK 

Poem JOC J)=NUJO(T) 

moo | +NUVEJOF ] 

NUPLJO( J) =NUJO(NUB) 

meMeEMB (J) =J 

KINDME(JU)=1 

SO TO "oY 

CONTINUE 

NO DIAGONAL(JOINT IS ON Pek hMe TER) 

i+ | 
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NUMIJ 


CONTINUE 

NUJ=BNUVE JU*NUMUIL =] 
NUJL=ENUVE JO+1 

K=)J 

DO 8 I=NUJ1sNUJ 
RATIO=VECOR(T)/VECOK (NUVEIUO) 
IF (RATIO=“-1le) 65270264 


CONTINUE 


THERE IS A DIAGONAL 
J=1-K+3*NUVEJUL*XNUHOUU-3*NUVL. 
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SUBRURT Ine 2H NAS 
DYNAMIC ARRAY HOCUR sVECUR s nNUIUs JURKENU RII U eWUPL UYU SNUG 1 UU ON UME 
OYNAMIC ARKAY K INGME »SUJO] sHOUCOCR Iie Tl 
COMMON HOLEN + VELEM., THICK »sMAMUIL » YOUNG s PRAT Tu sHOLEWN] s HULEWN2 3 VELEN] 
leVELENZsEXsWHYsK INP] EX 2 omHYe] sl-ULOM1 = SIT LILUsCl YP esi ieene eeee 
2NJI sAdUCUR(P) sVECUIK(P) sNUJU( PP? sYUKENU (FP?) srmiULLing 5s VEERING sINUR Uae 
Z3esNUVEIJ2ZeNUVEIJ3 eNUPLJU(P) eNUM[JU(P) NOME D (FE? eS INDMC (PY) sINOML eNO ome 
4 »sNUHOMB SNUXRIDI>s NULEUL-eNJUPULUL sHUdU] (FP) SHUCUSA LIP Ys VECOR] (ee 
PNULDS»ENDLD>LDMSITAsLOWMENDS IN Deore oe ee ey 
WOVEJGENUVEJLANUVEU 2 FNUV EJS 
NUHOJO=NUHOJS14+NUHOU24+NUHOJ 3-2 
Seer ORMAT (315) 
DO 4 [T=15NUME 
iNew: ( l)=u 
4 CONTINUE 
G VMERTICAL “EMBERS 
NUJENJI-1 
K=u 
DO 2 T=1sNUJ 
NAL PHA=NUVEJO#NUHOJ1+(NUVEJ1—1) 
Pee TA=NALPHA+ (NUROJ2=3)* (NUVEIIFRUV Eos 
IF (I-]-NALPHA) 12212511 
Pee CONTINUE 
Pees T—-NBETA) pipe en Wr be 
me CONTINVE 
a RBPEGULAR PROCEDURE 
meee DUCCR (1) —HOCOR( 141)) L225 02a 
fee CONTINVJE 
wie 1 —K 
eee J) =NUICt L) 
mereel JU J)eNUJOCIF)T} 
NUMEMBE (J) =J 
SO TO 2v 
fee CONTINUE 
K=K+1 
SO TO 2h 
me CONTINUE 
fee (VECGCOR( 1)-VECOR(NUVEJI1?7) Il3 lees 
me CONTINUE 
ee 10 i2 
mo CONTINUE 
& Port Ze«MEMBERS 
NALPHA=NALPHA-(NUVEJO-—1L) 
NUJ=NIJ-NUVEJO 
K=u 
mons | =e NU 
tr (1=NALCPHA) 22322321 
mee CONTINUE 
le (lHNBETAM@ 23322322 
mee CONTINUE 
C REGULAR PROCEDURE 
J=NUVEMS+ 1K 
NUM JOC) =SNUU0 CT) 
NSUb=I+4NUVEJO 





























Ze 


Zoi. 





Nel Low ov) SNe ( § Be 

mime Mot J) =2 

GeO 32 

CONTINUE 

IF (VECOR( IT )-VECUR(BUVEJ1)) 23537%545e 3% 
CONTINUE 

NUB=NUVEJ1L+NUVEJ2-] 

mee 6 VECOR (1 J—-VECOR(UNUB) ) 9252525 ees 
CONTINUE 

K=K+1 

oe TO 3 

MeN T INJE 

JENUVEMb+1-K 
NUP=I+NUVEJO-(NUVEJ2Z-2) 

NUMI JO( J) =NUJU(T) 

NUPLJU(J) =NUJO(NUP ) 

NUMEMB(J)=J 

Gen tO 3C 

BON TINUec 

RIGHT DIAGONAL “EMRERS 
NALPHA=NALPHA—-1 

NUJ=NUJ-1 

K=u 

DO 40 I[T=1,NUJ 

Meee I-NALPHA) 32:32:31 

@enw?;r I NUE 

Meee l—=NbDETA) 33932352 

CONTINUE 

REGULAR PROCEDURE 
NUP=I+NUVEJO+1 

MON TINUE 

mee OR (I )-VECOR(NUVEJO) J) 321.327 3s521 
CONTINUE 

J=NUVEMB+NUHOMB+ I-K 

Meet JO(J)=NUJO(I) 
NUPLJU( J) =NUJU(INUP } 

Rie rts( J) =J 

KINDME(J)=1 

See 1O 4u 

CONTINUE 

K=K+1] 

GO TO 4v 

fom) I NUE 

memeriNl TY OF HOLE 

Meet VECOR(1)-VtCOR(NUVECJ1)) 3352.33) 2551 
CONTINUE 

NUB=NUVEJI1+NUVEJ2-2 

Meee VECOR( I )—-VECUR(NUB) ) 32253229332 
Men? TINUE 

NUP=I4+NUVEJO-NUVEJ2+3 

fem tO 323 

Sen) INVUE 

LEFT DIAGONAL MEMBERS 

MeL PHA=NUVEJOFNUHO J14+CNUVEJL=1) 
NBETA=NBbETA+NUVEJO 
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PII = N= 1 
NUJTLENUVEJO+] 
K=u 


DO 5t I=NYJI1sNUJ 


TP WL NAWe EA) 42s 4D ste] 
CO Laie 

ie N si A 4 3 et Dowd? 
CONTINUE 


SeoUlLARK PROGEBEIRE 
Pie | —NU Vibe Ot) 


CONTINUE 

Mee (VECOR(I)-VECUR(NUVEJO) ) 421s Qa 
CONTINUE 
J=NUVEMB+NJHOML+NURIDI-NUVEJO4 I=K 

ieee 1 E (633) Locher 


Mee Ul J) =NUJO(T) 

Pete JU( J) =NUJOCNUP ) 
NUME™MB( JU) =J 

KINDME(J)=1 

GO, TO Su 

CONTINUE 

K=K+1 

mem 1O Se. 

@ernT I Nve 

Pee NITY OF HOLE 

Meee EGOix (1 )-VECOR(NUVEJI)) 
CONTINUE 
NUB=NUVEJ1L+NUVEJ2-2 

Meee VECOK (I )—-VLCurR( NUL) ) 
CONTINUE 
NUP=I-NUVEJOF+NUVEJ2-1 
GOs @ 423 

CONTINUE 

RETURN 

END 
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SUBROUTINE POBOC 
DYNAMIC ARRAY HOCOR»sVECOR sNUJOsJORENU sKINDJOsNUPL JO»NUMIJO»NUMEMB 
DYNAMIC ARRAY KINDMEsNUJO1sHOCORI »VECORI 
COMMON HOLEN>s VELEN» THICK »sMANUJO»s YOUNG sPRATIOsHOLENI] »sHOLEN2»sVELEN1. - 
L»sVELEN2:sEX sWHY sKINDI 29EX29WHY2 sPOLOM1 »sKILILO»sLTYPE »NUHOJI1 »sNUVEUJ]1 » 

BS  2NIJsHOCOR(P) »sVECOR(P) sNUJO(P) +s JORENULP ) sHOLFEN3 sVELEN3 sNUHOI? sNUHOI3 
3sNUVEJ2 sNUVEJ3 sNUPLJO(P) »sNUMIJO(P) sNUMEMB(P) sKINDME(P) »NUME sNUVEMB 

4 4s+NUHOMB »sNURIDI sNULEDI sNUPOLOs NUJO1] (P) sHOCORIT(P) sVECORI(P) sNSIZEs — 
SNULDSsENDLDsLDMSTAsLDMENDsKINDJO(P) sK IOUT »sNLOTY 


— — A gare eetatn ee = —_———— 





























Cc THIS_ SUBROUTINE TAKFS EACH PARTCULAR POINT SUPPORT <WiTH — 
C COORDINATESsEX( I) AND WHY( IT) AND OF KINDsKIND1I( I) AND MATCHES 
TH THE NE ARE S TG RD NE TT Fe RET Fea ee 
c NUJOCALPHA) WITH COORDe HOCOR(ALPHA) AND VECOR( ALPHA) 
Peo itsetFX —HOCOR 1) ) Aer ( WHY. =A ECORI) ) lca i 
NALPHA=]1 
= om lOO JEP»N) ee a _ en 
BIST=(EX -HOCOR (J) )**2+( WHY —-VECOR(J))**2 


DiS T_1 — Defoe) Oe 1 O 0 30 O98 ee 
Le CONT INUE 
—__DILSTANCE FROM NEW POINT GREATER THAN OR EQUAL TO PREVIOUS 
Go 1O 100 
ee ONT INUE_ ee 
a DISTANCE FROM NEW POINT LESS THAN PREVIOUS 
3) =) 
MAEPHA=J 
mee TOOSCONTINUF 
aq NEAREST JOINT FOUNT 
IF (KINDJO(NALPHA)-KIND1 fli oO Ors 2:00 
“150 KINDJO(NALPHA)=KINDI1 
DOLGON TINUE 
RETURN 
= END, a  __. ne 
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SUBROUTINE LIBOC 


111, 


DYNAMIC ARRAY HOCORsVECOR sNUJO: JORENUs KINDJO»sNUPL JO sNUMI JO sNUMEMB 


DYNAMIC ARRAY KINDMEsNUJO1»HOCOR1 »sVECORIL 


COMMON HOLEN»sVELEN » THICK» MANUJO » YOUNG » PRATIO»sHOLEN1»>HOLEN22VELEN1 
1» VELEN2sEX sWHY sKIND1 9 EX29 WHY 2 sPOLOM1] sKILILO»LTYPE sNUHOJ1 sNUVEJ1> 
2NJ »HOCOR(P) sVECOR(P) »NUJO(P) » JORENU(P) s HOLEN3 » VELEN3 sNUHOJ2 » NUHOJ3 


ae 
eS 


3 sNUVEJ2 sNUVEJ3 sNUPLUO(P) sNUMIJO(P) sNUMEMB(P) sKINDME(P) »NUME »>NUVEMB 
4 » NUHOMB »NURIDI »NULEDI »NUPOLOsNUJO1(P) sHOCOR1(P) sVECORI1(P) »NSIZE> 





SNULDSsENDLDsLDMSTAsLDMENDsKINDJOCP) »sKIOUT sNLOTY 
EQUIVALENCE (EX2EX1)2(WHY»WHY1) s(KIND1sKIND2) 





“KG THIS SUBROUTINE TAKES A SET OF LINE SUPPORTS AND APPROXIMATES 
_¢ THEM BY A SET OF EQUIVALENT GRID JOINT SUPPORS 


PO 100 J=lsNJ 
_¢ FIND EXACTLY MATCHING COOR. 
IF (EX ~EX2 ) 20510320 
4 10 CONTINUE __ _ — 
a X=COORDe MATCHES 
IF (EX1 ~HOCOR(J)) 11912911 
11 CONTINUE 
@ NOT A JOINT ON BOUNDARY | 
GO TO 100 
12 CONTINUE 




















€ A JOINT ON BOUNDARY 
IF (WHY2 — -VECOR(J))13914914 
13 CONTINUE 
iro JOINT OUTSIDE Y-RANGE 
GOr10" 100 
— -14 CONTINUE 
C JOINT PROBABLY INSIDE Y RANGE sCHECK Yle | 


| IF (WHY1 -VECOR(J)) 16216915 
15 CONTINUE 
: 























C «JOINT OUTSIDE Y-RANGE 
( Soro Too 
«16 CONTINUE 
ec JOINT DEFINITELY WITHIN Y-RANGE 
JORENU( J) =JORENU(J) +1 
GO TO 80 
20 CONTINUE — 
a Y-COORDe MATCHES 
| IF (WHY1  -VECOR(J)) 21522921 _ 
21 CONTINUE 
C NOT A JOINT ON BOUNDARY 
GO TO 100 
(22 CONTINUE | 
C A JOINT ON BOUNDARY 
EX HOCOR (b) 23 9 2duuod 
23 CONTINUE. 
E JOINT OUTSIDE X-RANGE 
SO TO 100 
24 CONTINUE. —— 
e JOINT PROBABLY INSIDE X-RANGE»CHECK X1 
ie (6%)  =HdOeOR(iye2g82E soe 























| 25 CONTINUE 
C JOINT OUTSIDE X-RANGE 
GO TO 100 
26 CONTINUE . ——e 
JORENU(J) =JORENU( J) +2 
C «JOINT DEFINITELY WITHIN X-RANGE 







































° z S 
«ot i ' ( 7 1) \ \ ° ¥ 
Write iV ‘ - z - i \ 5 ; ¥ - 
ae. 6 A ‘1 os if . ; . —' rrr, Age 
oa} “Wes ; » ote) = Feeait . 4 ee ee ee ee 
mae | " ; ‘ , l» “ i — ti eees 
O40 heta “ae ale i-e@ 65-090). vee sia laa 
Py I ; Die itl ss 1a) @ ies Tie — & yey 
— sie “te i=) oo © rn’ ; 94 ee 
*..¢ _— 
he ae’ tte Ras 
(te @ : "ai- : 1 
| ntl 
Po -J 
oAji- ' AAs YT vat ant 
et le ra 
er 
| 2 oo bon rie 
ut ela eee ral 
res by 
Sh te on het BY 
>i 
we a 
- ov): a \ eles (co 
re ae 
_ i 


ee 
— 


aly Waite Farle 
1?) s bees Mee | 


ar 
(Soave ae Lib tee et- 


as a at 
hat hap er Se 


Frere lll a tT 
ui | oth 
wai if 

- » F She | Laut? & magia ng 
Stic 10001 he ‘at 
mF ae 
“ad Bam @ Lary ™* y. 
18? (oe 


tio: ee ae ey. ie oe 
wa ly). + pipe TaD 





ee 


eee SO TF (KINDS = Os10051]00. 
90 WI NDJO (ok I ND? 
100 CONTINUE 


















































C PICK ANOTHER JOINT = 7 
200 CONTINUE ————_— 
C PICK ANOTHER BOUNDARY CONDITION 
PRR TUM 7 
END 
4 ae = a 
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mm ** 0 26. ——— 
SUBROUTINE POINLO 
3 DYNAMIC. ARRAY HOCOR sVECOR sNUJOs,sJORENU sKINDIUOsNUPL JOsNUMIJOsNUMEMB 
DYNAMIC ARRAY KINDMEsNUJO1] sHOCOR1»sVECORI 
a COMMON HOLENs VELEN» THICK »MANUJO s YOUNG» PRAT 1OsHOLEN1»sHOLEN2» VELEN] 
lsVELEN2sEX sWHY sKIND1 sEX29WHY2 sPOLOM] sKILILOsLTYPE sNUHOJ1 sNUVEJI1 > 
—2NJsHOCOR(P) »sVECOR(P) sNUJO(P) »JORENU(P ) sHOLFN3s VELEN aNUHOIJ2 sNUHOJ3 _ 
3 »sNUVEJ2 sNUVEJ3 »sNUPLJO(P) »NUMI JO(P) sNUMEMB(OP) sKINDME(P) sNUME s5NUVEMB 
4 4sNUHOMBsNURIDI sNULEDI sNUPOLOsNUJO1 (P) sHOCORI(P) »VECORI(P) »sNSIZE>5 —_— 
SNULDSsENDLDsLDMSTAsLDMEND sKINDJO(P) sK IOUT sNLOTY 
_ FQUIVALENCE JEXsEXLOAD*»)WHYsWYLOAD* 
a fais OUBROUTINE TAKES EACH POINT sacecnm AND FINDS THE NEAREST NODE 
feo THF GRID WHICH JS NOT A SUPPORTsTO BE CONCIDERED AS POO ee 
@ APPLICATION OF LOAD. 
eee OR MAT OT 86) 2H JOINTL-tOADS) a ee = 
feroRMAT (T83I145s9H FORCE Z sF15e3) 
momPORMAT (T8si4s10H MOMENT X 5F15,3) 
PemeORMAL (T8sl1H LIST LOADS 314} 
: ORMA R OH MOMEN a3) 
SemrORMAT (T8s20H LIST DISPLACEMENTS 514) 
___ FORM AN ARRAY OF NON-SUPPORT JOINTS. _ >... 
CaLk NOSUJO 
eee (NUPOILO—1). 22521922 
Pe CONTINUE 
NK efQ) EEE Ee 
22 eee ee ey OO ee eee 
LPHA=1 em 
ae foo  J= INST CE 
fe = -EXLOAD-HOCOR] (J) ) **2+( WYLOAD=VECOR] (J) ——F7? EEE 
Mm tDISTI~DIST) 1003100330 
OinD =12 
NALPHA=J 
m= 100 CONTINUE. _— 
a MeAREST JOINT FOUND 
ere (10) 5101510351035 1032106)sKIOUT 
ioremeONT TINUE 
PRPOS=NUIJOTI(NALPHA) a 
WRITE (7383) NUJOS 
eto TO 200 _ 2 — 
106 CONTINUE 
pee NU JO9S=NUJO] (NAL PHA). _ = 
Merrie (7386) NUJOS9 
e 0 OU 
101 CONTINUE 
mes | O.LY~)) 9229} 297 ne ee eee 
91 CONTINUE 
ema CONTAINS AN }el OAD JS A Z FORCE 7s —orr ’ —_ Tae 
NUJO9S=NUJO1 (NALPHA ) 
R L7 wf 5) NUJOSOsPOLOM) aaa aaa 
=o 10 200 
mee? JC ONT I NUE : ee ee ee KK 
re NLOTY DOESN'T CONTAIN AN JeCHECK IF 2 DOR ie 
ee ONE OLY=2). 94993294 ~ «x - - — ——_—_—_ 
93 CONTINUE 
mt EY CONTAINS A 2g 
NUJO9S=NUJO1 (NALPHA) 
- Pe Pe (eB O) NUJOOsPOLOM)] >? eee 4 
GO TO 200 
meee CONT LNUE WW === Se > 
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NUJO9=NUJO1 (NALPHA) 
- WRITE (7985) NUJO9 sPOLOM1 — ill 
200 CONTINUE 
So LCMUNUPOLO=NUPOLO#L — 
RETURN 
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= 2% 0 2 6. 2 2a 
SUBROUTINE NOSUJO 
DYNAMIC ARRAY. HOCORsVECORsNUJOs JORENU sKINDJO»sNUPLJO»s NUMI JO »NUMEMB 
DYNAMIC ARRAY KINDMEsNUJO1sHOCOR] »VECOR]I 
COMMON HOLENs VELEN s THICK »MANUJOs YOUNGs PRATIOsHOLEN1] sHOLEN2sVELENL 
LlesVELEN2 sEXeWHY sKIND1 2s EX29WHY2 sPOLOM] sKILILOsLTYPE sNUHOJ1 »sNUVEJ]1 » 
—___...._. 2NJ»sHOCOR(P) sVECOR(P) sNUJO(P ) » JORENU(P ) »>HOLFN3sVELEN3 eNUHOJ2 sNUHOJ3 
3 sxNUVEJ2 »sNUVEJ3 sNUPLJO(P) sNUMIJO(P) sNUMEMB(P) »sKINDME(P) » NUME sNUVEMB 
—__ 4sNUHOMB»NURIDI sNULEDI sNUPOLO»sNUJO1(P) sHOCORI(P) sVECORI(P) sNSIZE> 
SNULDS sENDLED»sLDMSTA sLDMENDsKINDJO(P) sKIOUT»sNLOTY 
ee -KSO 
me 50 TE TsNJ 
ee et CK INDIUO(CT ) =O) TOsSeq70 LL EEE 
D PCONTINUE 
_— NOT..A_ SUPPORT _ a 
J=I-K 
ee Ll) ) SNL IO CT) a ee 
MOCORI(J)=HOCOR(] ) 
Bee OR) ( J)=EVECOR(T) OO —=E 
ve To SO 
MOSEONTINUE — as See 
a meer PORT JOINT SKIP. 
Saati i. = 
50 CONTINUE 
ee a 
RETURN 
END 
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BeXCO2 Gee 
*¥*OBJECT 
SUBROUTINE. JOINRE : 
DYNAMIC ARRAY HOCOR sVECOR »NUJO»s JORENU sKINDJOs NUPLJOsNUMI JO sNUMEMB 
= _ DYNAMIC. ARRAY KINDME»sNUJO1sHOCOR1sVECOR] ~ = 
DYNAMIC ARRAY NU 
_________ COMMON. HOLENsVELENs THICK sMANUJOsYOUNGsPRATIOsHOLEN]sHOLEN2,VEILEN] 
ls VELEN2 sEX sWHY sKIND1 9 EX2sWHY 2 sPOLOM] sKILILOsLTYPE sNUHOJ]1 sNUVEJ1 » 
q 2NJ sHOCOR(P) sVECOR(P) sNUJO(P) » JORENU(P) »>HOLEN3 »VELEN3 sNUHOJ2 sNUHOJ3 
3 sNUVEJ2 sNUVEJ3 sNUPLJO(P) sNUMIJO(P) sNUMEMB(P) sKINDME(P) sNUME sNUVEMB 
3 4»NUHOMB »>NURIDI sNULEDI »NUPOLO»sNUJO1(P) sHOCORI(P) sVECOR1(P) »sNSIZE > 
S5NULDSsENDLDsLDMSTAsLDMENDsKINDJO(P) sKIOUTsNLOTY 
Sp 3 FORMAT (T8,20H TYPE PLANF GRID xy ) 
4 FORMAT (T8s15H CONSTANTS E 9» sFlle3s5HsALL ) 
a SS eORMAT oT 8e15H CONSTANTS G 2 sFlle3s5HeALL J. 
6 FORMAT (T8s18HJOINT COORDINATES ) 
8 FORMAT (7T891422F10e4s9H SUPPORT ) 
10 FORMAT (T821492F10e4) 
| 12 FORMAT (T8s15HJOINT RELEASES ) 
14 FORMAT (T8s1I14912H MOMENT X Y ) 
15 FORMAT (T8sI4210H MOMENT X_) 
16 FORMAT (T8514s10H MOMENT Y ) 
___17 FORMAT (T819H UNITS» INCHES »sKIPS ? 
WRITE (7:93) 
Meeteetye4) YOUNG 
G=YOUNG/(26*(1e+PRATIO) ) 
— WRITE (725) G SS 
WRITE (7917) 
WRITE ( 7is6) 
oe] 
SS ON Oe 
Be 100 [=lsNJ 
KI EK INDJOCT) +1 a 
GO TO (80350370550) sKI 
eee GON INUE ae, 
q JOINT IS HINGED»sRELEASE MOMENTS 
OU NE 
eal 
NUJI=NUJO(T ) 
HOCO9=HOCOR(I ) 
_ VECO9=VECOR(I) 
WRITE (7598) NUJ9sHOCO9sVECO9 
eoe70.300 
70 CONTIPUE 
<< Seen 1 SF TED N.S 
NUJ9=NUJO(T) 
meOO=HOCOR( J) ie 
VECO9=VECOR(T) 
WRITE (758) NUJ9sHOCO9 »sVECO9 
GO TO 100 
80 CONTINUE a go a EEE 
' Ie JOINT IS FREE 
~NUJI=ENUJO(T) 
HOCO9=HOCOR(T) 
VECO9=VECOR (I) 
WRITE (7910) NUJ92HOCO9 sVECO9 
PTOO CONTINUE  __ 'd . —— eee hl! 
. IF (K=-1) 110591105120 
mm 110 CONTINUE | _ _ —— 
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THERE ARE NO. JOTNi wR E bee ASE Se eee 
GO” TOL U0 

CONTINUE . 

THERE ARE K=-1 JOINT RELEASES 
Vir ~ (ee) 

L=K-] 

BO—-) SO" K = 15h ae 
NOCO=QQFDIL(NUs«K ) 
KI=QQFDI1(JORENUsNOCO) +1 a. i 
oe 10 (1470993059730) skK1 








ee gn! SOLCONT INUE P : “ 


C 


A LINE SUPPORTsRELEASE ONLY ONF MOMENT 


EE SI © 


132 


———— 


140 


m0 


JOR=JORENU(NALPHA) 

feet (131 S32 ee JOR OO O— te a - - _ = 
NO=NU(K) 

WRITE (7516) NO. =< = 

oe TO 150 

ee ee 
WRITE (75915) NO 








mO510 150 a 


NO=NU(K ) 
WRITE (7914) NO 





CONTINUE 


eo eeNT LNUE 


pn | 


202 


a ee 


300 CONTINUE 
WPPSTROY. NU. ———— eT 


FORMAT (T8s19H MEMBER INCIDENCES ) 

Merl (75201) _—_ ee ee — as — —— 
FORMAT (T8s4I14) 
DO 300 I=1lsNUME 














NUMEMS=NUMEMB(T) 
NUPLJIJO=NUPLJUO(]) 
NUMI J9=NUMIJO(T) 
WRITE (75202) NUMEMS »NUMIJ9»sNUPLJO 








RETURN 
ND Ss 
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ete 6 ny a ee eee ee 
SUBROUTINE MEMPRO 
DYNAMIC ARRAY HOCOR>sVECORsNUJOsJORENUsKINDJOsNUPL UO »NUMI JO »NUMEMB 
DYNAMIC ARRAY KINDMFE »sNUJO1sHOCORL»sVECORI 
DYNAMIC ARRAY MEMNUsREINERsTOINER | 
COMMON HOLENsVELENsTHICKs MANUJOs YOUNG sPRATIOsHOLEN1 sHOLEN2 s VELEN 
———-—J-»VELEN2sEXsWHY sKIND] »EX2 »sWHY 2 sPOLOMI] sKItLIlLOst TYPE sNUHOIUL »s NUVE J}. _- 
2NJ sHOCOR(P) sVECOR(P) sNUJO(P) s JORENU(P) sHOLEN3 9 VELEN3 sNUHOJ2 s NUHOU3 
7. 3»NUVEJ2 sNUVEJ3 sNUPLJO(P) sNUMI JO(P) sNUMEMB(P) sKINDME(P) sNUMEsNUVEMB  __- 
49 NUHOMB »sNURIDI »sNULEDI »sNUPOLO»sNUJO1 (P) sHOCOR1I(P) sVECORI(P) sNSIZE > 
ed ONULED Ss» ENDLD > 1 DMSTA »L_DMENDs KINDUGEP eK TOUT »NLOTY 
DEFINE MEMNU »NUME 
ee DEFINE REINEReNUME 
DEFINE TOINER»sNUME 
ee FORMAT (T8s28HMEMBER PROPFRTIIESsPRESMATTC )  |o—=§ = 2. 
Se evAT (T8elte4H LY »F1l56¢S9%4H 1X, 9 S59 
eo FORMAT. (7180213H LOADING "2" )_ ae 
FACTOR=(THICK¥*3) /(24e*(1le—-PRATIO¥*#2) ) 
Peeereitc(f-e) Doe ee  CCCCC“‘(§HUUUUOCté‘(‘#WWWtC(‘( KUO ica 
B® 120 [T=] +NUME 




















- HE LNERA I) =O. 
10 TOINER( 1) =0 
mc...___ ARRAYS SET FQUAL TO ZERO 











I 1 =NUHOMB+NUVEMB 
=T1+NURIDI 
Be 200 T=I12I2 
pe O 
be yoo J= sil 1. 
NALPHA=NUMIJOCT) 
NBETA=NUPLJO(T ) 
NA PH=NUMI JO ( wr) 
NBET=NUPLJO(J) 
 _ TF (NALPHA=NALPH) 1009503100 
50 N=N+] 
__ MEMNU(N) =NUMEMB(J) 
100 CONTINUE 
= pa 
ee 'O (20032002110) sNAO 
110 CONTINUE 
emi cO K=lsI1 
_NALPH=NUMIJO(K) 
NBET=NUPLJO(CK ) 
(NBETA-NBET) 1201113120 
met N=N+)] 
MEMNU(N)=NUMEMB(K) . _ 
fee ONT INUE 
Meee SO Kael2eoNUME a 
NALPH=NUMIJO(K) 
NBET=NUPLJO(K ) 
NCAT=MEMNU(] ) 
IF (NBET-NUPL JO(NCAT)) 13091219130 










































121 N=N+] 
____ MEMNU(N) =NUMEMBCK) i : — ae 
130 en NUE 
=N+1 = 
. MEMNUTN) SNUMEMBCT) 
ie SIX BEAMS OF PANNEL FOUND 





VESIDE=VECOR(NBETA) -VECOR(NALPHA) 
HOSI DE=HOCOR (NBETA) -HOCOR(NALPHA) — 
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fe ADT AGSSORTIVE ST DE REO+HOS IDE R* 2) EEE eee 
HOVERA=HOSIDE/VESIDE 
DIVERA=ADIAG/VESIDE 
DO'S Ouck = 15h) 
NGAMMA=MEMNU (CK) 
GO SO 161 5 16261 6 led 6251635063) 3K 
me) COR FIR HOVERA-PRAT 1O/HOV F Rife) 7 eS TG 
COEF 2=(2*HOVERA*(1¢+PRATIO) *(1.6-3e*PRATIO) )¥VESIDE 
6 O..).O- 17.0: — — “ So ee 
162 COEFI1=(1e—-(HOVERA**2) *PRATIO) *¥VESIDE 
——e COEF 2Z=(20*(1-+PRATIO)¥(1-¢-3e*PRATIO))*VESIDE iy % 
GO TO 170 
M™ 163 C = % ¥%% * 
COEF 2=0 
170 REINER(NGAMMA) =COEFI*FACTOR+REINER(NGAMMA) ss 
TOINER(NGAMMA )=COEF2*FACTOR+TOINER(NGAMMA) 
Bee GT INUE ee cel 
200 CONTINUE 
bie PON KegeaepCM  E—E—E aaa eee 
NUMEM9=NUMEMB(K) 
meme SQ=REF INER( K ) “ge 
TOINEQ=TOINER(K) 
NRITE (726) NUMEMS sRETINE9sTOINES —————————————————————— 
300 CONTINUE 
b RO ME MAN 
ODesetROY REINER 
DmolROY TOTNER 
WRITE (798) 
Se SREMRN eo ee 
END 
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me  %XO276 ws 
SUBROUTINE LINELO 
DYNAMIC ARRAY HOCOR »VECOR sNUJO»s JORENUsKINDJOsNUPL JOsNUMI JO »>NUMEMB. 
DYNAMIC ARRAY KINDMEsNUJO1sHOCOR1+sVECOR1 
_ _— DYNAMIC ARRAY XLOCAL »YLOCAL »sXGLLDsYGLLDsNULOMF >» KINDLOsNULOJO : ay 
COMMON HOLENsVELENs THICK sMANUJUO s YOUNG »s PRAT IOsHOLEN] sHOLEN2 sVELEN1 
oe UCL VEL EN2@ EX o WHY 0K INDI Ss EXD eo WHY? »POLOM] sK It Ti Oe! TYPE sO. |] 5 Nis ese 
Z2NJ sSHOCOR(P) sVECOR(P) sNUJO(P) s JORENU(P ) sHOLEN39VELEN3 sNUHOJ2 » NUHOU3 
a 3sNUVEJ2 sNUVEJ3 sNUPLJO(P) »>NUMIJO(P) sNUMEMB(P) sKINDME(P) »>NUME sNUVEMB 
4»NUHOMBsNURIDI sNULEDI sNUPOLOsNUJO1 (PP) sHOCORI(P) sVECORI(P) sNSIZE> 
Ss  ONULDS»sENDLD»t DMSTA»sLDMEND sKINDJO(P) sKIOULsNLOTY _. ; 
EQUIVALENCE (EXsXA) »(EX25XB) 9s (WHY sYA) s(WHY29YB) 9 (STALD»sPOLOM}1 ) 
eee OPE ENF XOCAL §sNUVEMBsNUHOMB a re ae ere 
DEFINE YLOCAL sNUVEMB+NUHOMB 
a DEFINE XGLLD  »NUVEMB+NUHOMB 
DEFINE YGLLD »sNUVEMB+NUHOMB 
C(O DEF INE NULOME sNUVEMB+NUHOMB ——————— 


DEFINE NULOJO »sNUVEMB+NUHOMB 
DEFINE KINDLOseNUVEMB+NUHOMB 


j K=0 
(a ——— 
Pee OXA=XB) 40541540 
__ 40 CONTINUE | 
IF (YA=-YB) 49242949 
ON EENVE 
G MeADINGSLINE PARALLEL TO Y- AXIS 
ee ees eS LN) ON *t fea 
IP (XA—-HOCOR(I)) 41154125411 
6) 1 CONT NUE. 7" _ ” —— aang 
C BOADING LINE DOESN'T COINCIDE WITH A SET OF MEMBERS 
K =(0 
N=0 
= G@» TO.100 ——— Be Aa 4 
412 CONTINUE 
ee merDING J INE COINCIDES WITH JA SET OF MEMBERS 
K=0 
Sere, JJ=)eNUVEMR ci 
NALPHA=NUMIJO(J) 
Seer rA=NUPLJO( J) 
IF (HOCOR(NALPHA)=-XA) 41994139419 
wees GONT INUE. : = <= —————rt 
IF (VECOR(NBETA)<-YA) 41994199414 
IN \ 
IP ¢VECOR(NALPHA) YB) 41524193419 
Ae, GOON Ge es Ba = 
IF (VECOR(NALPHA)-YA) 41654165417 
ene T LNUE —_ Se ee —— 
K=K+] 
XLOCAL(K) =YA-VECOR (NALPHA ) 
XGLLD(K)=YA 
=—@o FO 476 a 
417 CONTINUE 
K=K+1 
RIOCAL(KY=0 
XGLLD(K)=VECOR (NALPHA ) 
418 CONTINUE 
| NULOME(K)=NUMEMB(J) | , = 
KINDLO(K) =2 
ga IF (VECOR(NBETA)~YB) 4161°9416194162 
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S436) CONTINUE ee a 
YGLLD(K)=VECOR(NBETA) 
YLOCAL(K) =VECOR(NBETA)-VECOR(NALPHA) am ~ a oe 
GO TO 419 

— 4162 CONTINUE 
YGLLD(K)=YB 

ee VOCAL LK ) = Y Buel F CORN NDA) aS eee 

419 CONTINUE 
> GO _TO. 200 | — — ES 
42 CONTINUE 

- C LOADING. LINE. PARALLEL TO X-AXIS a 
DO 420) Tt =aesiny 

es ee ( Yo VE COR (1) ) 427) ote22—e49) eee 

421 CONTINUE 
GOxTO49 _ ——e 
422 CONTINUE 

feo TOAD INY INF COINCIDES WITH A SET OF MEMBERS 

| NUV=NUVEMB4] 





























DO 429 J=NUVsNUVE 
NALPHA=NUMIJO( J) 
NBETA=NUPLJO( J) 

















— im ( VECOR(NAL PHA) =-YA) 429:2423.,429 
423 CONTINUE 
MOCOR (NE A)J—-XA 2 5 





424 CONTINUE 

Meme(HOCOR (NAL PHA)=XB) 42594295429 
ez aeNTT NUE 

Pee bHOCOR ONAMPHAY —XArG2 624269427 
426 CONTINUE 


a, eee EN 


XLOCAL(K)=XA—-HOCOR ( NALPHA ) 
X Deere XA 
GO TO 428 
4 Meee NWI 
K=K+] 
meOoCc AIK) =0 
MGWED (K) =HOCOR (NALPHA) 
MeemoemmINUE oC 
NULOME (K)=NUMEMB( J) 
aC KTNDLO( KDE? ; one aoe 
eee (HOCORMNBETA)—-XB) 4261°2426134262 
Mee i 
meee D(K )=HOCOR (NBETA) 
meeOGrts (K )SHOCOR(NBETA)—-HOCOR(NALPHA) oo 
GO TO 429 
62 GONTINUE | ————————————— Se ee 
YGLLD(K)=XB 
Oley: K )=XB-HOCOR (NAL PHA 
429 CONTINUE 
OTL. 200. inn, MN TIN —— 
49 CONTINUE 
oe FP OADING | INE NOT_OF_THE FORM X=CONSTE sOR Y=CONSTg= — 
N=0 
—_ ae 
COEF=(YB~-YA)/(XB-XA) 
=e DO DO D=.) sw NW SLE MBs 
2 NALPHA=NUMIJOCT) 
7 AE le eNO ee Oe ee = 4c nn ae... 
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aL FON = 04) fyi 9S ph Ee 
494 CONTINUE 
= DO 493 J=15N 
IF (NALPHA=NULOJO(N) ) 4935100,493 
a 429.3..CON Tel NUE 
Y=YA+COFF *(HOCOR( NALPHA) -—XA) 
= ——_—______ _TFEDP!|=ABS(VECOR(NBETA)—y) 
TED=(VECOR(NBETA) -VECOR( NALPHA))/1000. 
ABSL=ABS LY)= : = 
IF (VECOR(NALPHA)-0O) 49134965491 
= Fae fle al NIE ‘ _ _ = 
me A ABSY=1TEO 5755753491 
eumemmeeemmeee NGS Sa Eee eee 
IF (TEDPL=-TFD) 563562492 
oe 4? CONTINUE. : —— 
RPmPOKARKHOCOR(NALPHA)) 509505100 
| SO -GON TI NUE 



































ie (XB= HOGOR(NALPHA)) 100651551 
ee GONT UN 
C INTERSECTION POSSIBLE BUT NOT CERTAIN 


IF SVECOR(NAL PHA) -Y) 5221002100 | 
52 GONTINVE 
eee re EC TLON, STIL POSSIBLE 
Bee eve COR CUNBE TAO™=Y) 100556555 
ONT LNUWE 
c MHTERSECTION CERTAIN 
K=Kt] i, es 
NULOME(K)=NUMEMB(T) 
~i> KMEOGAL LK ZY-VECOR(NALPHA) 
XGLLD(K)=HOCOR (NALPHA ) 
0) 
KINDLO(K) =0 
"GGn TO 1.00 
56 CONTINUE 
eee JO WNT »~OF MEMBER &LhES ON EGAD LIME 
PP CK INDYOONBETA) -)) 56)e5629562 
OS Gr i ee 
© JOINT NOT A SUPPORT 
_ Kektl] _ sea 3 
NULOME(K)=NBETA 
_ XGLLD(K)=HOCOR(NBETA) 
YGLLD(K)=VECOR(NBETA) 
Bene) mr a 
N=N+1 
NULOJO(N)=NBETA 
IF (YA-YB) 10032005100 
Poo? Ne] __ 2 a ee 
NULOJO(N) =NBETA 
— Ge; TO _.1.00 
57 CONTINUE 
























































eo MINUS JOINT PIES ON LOAD LINE 
IF (KINDJO(NALPHA)=-1) 57195729572 
571 CONTINUE a... ——s — 
a JOINT NOT A SUPPORT 
K=K+] ———— 
NULOME (K) =NALPHA 
mee XGLLD (KY =HOCOR(NALPHA ). 





YGLLD(K) =VECOR (NALPHA) 
KINDLO(K)=1 _ 
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= a AE 
NULOJO(N)=NALPHA 
Lele AS eRe OnesiO'Dee,) OF) 
Dale GONDLNVUE 
4a N=N+1 
NULOJO(N) =NALPHA 
me UL O00 hUCONT INE 





C ALL INTERSECTION POINTS WITH VERTICALMEMBERS FOUNTeTHERE ARE ALSO 
=< AiO) New TER SEC Toh@nisi 7 
S HORIZONTAL MEMBERS 
a DpO-2 00ND EX= bw HOMB. 
I=INDEX+NUVEMB 
ee Ne P HAS NIIMT IOC Te eee eee eee 
NBETA=NUPLJO(I) 
Be eN=0 jee, lage —=Ee 
130 CONTINUE 
Reese lIN OF <THE | OADING, PaVSSF S) THRUSSO@ME JODNTS. ae 
DO 13°77 J=15N 
BA PANU Geb le. 5 lube Daiwa O 9 135 eee 
o> JGONT bNUE 
— eee NBETA-NULOJOCS) ) 13975200, 137 - 
ae CONT TINUE 
eee MEMBER J. DOESNOT PASS THRU A LOADED JOINT 
Mae 1 A~YB bul 38% 2005139 
- ONT IN 
IF (YB-~VECOR(NALPHA)) 2003138151381 
Ba GONDINUE  _ a ee 
Me CT YA-VECORCUNALPHA)) 151,751 5200 
eee NIE. .« — © wee 
IF (YA-VECOR(NALPHA)) 2003139131391 


SE ae oae see — ee ee 





























: 6 a 
BeOra-VECOR(NALPHA)) 1519151:;200 
EGOMMIENUF 





@ PAMeRSECTION POSSIBLE BUT NOT CERTAIN 
— X=XA+((VECOR(NALPHA)~YA)*(XB-XA) )/A(YB-YA 
IF (HOCOR(NALPHA)=-X) 15292003200 
ON N 
‘@ POMPE RSECTION+STIl&t POSSIBME 
Pere HOCORANBE TA)-X9 20020051755 OO “ee 
155 CONTINUE 
Ree PERSE CTION CERTADN: 
K=K+] 
N OME (K ) =NUMEMB( I ) 
YGLLD(K)=VECOR (NALPHA) 











pelea ( K)= 
XLOCAL(K)=X-HOCOR (NALPHA) 
EIGEN DEOWK ) =O ee 





200 CONTINUE 





\ = 
CALL LINE (S$KINDLOsSXGLLDsSYGLLDsS5NULOME s SXLOCAL s SYLOCAL ) 

— Bee uO Y meleOC Ain, ; —_—= — — — a 
Peo ROY YLOCAL 

foe 6 DPESTROY. XGLLD = - o- ———— 
mesoTROY YGLLD 
® adl@ OM 
DESTROY KINDLO 

_— pe STROLL OO : i er 
RETURN 

ee EF ND. —— a al —— aa 
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SUBROUTINE LINE(KINDLOsXGLLDsYGLLD»sNULOME »XLOCAL » YLOCAL) 

a DYNAMIC. ARRAY HOCOR sVECOR sNUJOsJORENU sKINDUOsNUPLUJOsNUMI JO, NUMEMB 
DYNAMIC ARRAY KINDMEsNUJO1sHOCOR1+VECORI 

a DYNAMIC ARRAY XLOCAL »sYLOCAL»sXGLLDsYGLLDsNULOME sKINDLO = 
DYNAMIC ARRAY RADIUSsSTALME oP sJOINT 

mC OMMON HOLLEN VELENeo THICK sMANUJOSYOUNGsPRATIOsHOlEN]sHOI EN? «VFICE Ns 
lsVELEN2 EX WHY sKIND1 2 EX2 sWHY2 sPOLOM1 sKILILOsLTYPE sNUHOJ1 sNUVEVJ] » 

_ 2NJsHOCOR(P) sVECOR(P) sNUJO(P) s JORENU(P) »HOLEN3 » VELEN3 sNUHCJ2 » NUHOJ3 
3»NUVEJ2»sNUVEJ3 sNUPLJO(P) sNUMI JO(P) sNUMEMB(P) sKINDMECP) sNUMEsNUVEMB 

—_—. 4sNUHOMB»sNURTDIsNULEDI sNUPOLOsNUJO1(P) sHOCORI(P) sVECORI1(P)»sNSIZEs5 
SNULDS sENDLD»sLDMSTAsLDMENDsKINDJO(P) »K IOUT »NLOTY 

ee MOLI TV AL ENCE (FXeXA)»o(EX2sXB) s (WHY sYA) » (WHY? sYB) s(STAl Ds POROM Ia 





























a THIS SUBROUTINE IS A LOGICAL CONTINUATION OF LINELOeIT WAS FORMED 
ete eee ASA SFRARATE PROGRAM SINCE PRECOMPILFER AT THIS TIME CANNOT ACCEPT 
a SUCH A LARGE NUMBER OF SYMBOLS 
 am- seORT “21ST OF POINTS ACCORDING TO THEIR DISTANCE FROM THE > 
ie PED PNG GDGE OF THE LINE LOAD 


Seer INE RADIUS s«NUVEMR+NUHOMBR cies see: i 
DEFINE STALME sNUVEMB+NUHOMB 
fee Pe] NEF Pe sNUVEMB+NUHOMB  _ —— ae 
PeriNE JO DNT sNUVEMB+NUHOMB 
fee 10 FORMAT (18s12HUWOINT LOADS) 7 Z 
PY MORMAT™(T8,h4,9H FORCE Z 3F15e3) 
ORMA Be Tl4910H MOMENT X 3F 1503) 
iS GORMAT (T8914s10H MOMENT Y oF 150e3) 
4 FORMAT(T8913H MEMBER LOADS) “ 
MomeermAl (18s143s23H FORCE Z CONCENTRATED PoF15e332H LsFl5e3) 
16 FORMAT (T8sI14924H MOMENT X CONCENTRATED PosF15e392H LoF1l5e3) 
17 FORMAT (T8s14s24H MOMENT Y CONCENTRATED P9F15e392H LoF15e3) 
poses _B SF15e3) 
SomeereMmAT (Theta .20H MOMENT _X UNIFORM oW oF] 5¢ Seve LA sFibs3aeH = 7 
1T894H LB 9F15e3) 
ee mORMAT (T38el4s20H MOMENT Y¥ 
1P6S4H LB oF P56¢3) 
ememMAT (T8e ]4e19H FORCE Z loINE AR WAmseil Sie Siete B coe EoleSie Aisgiee/ 
ippeeeeH WAY »F15-3s4GH LB »F150e3) 
22 FORMAT (1T85:14s20H MOMENT X LINEAR WA »9F15e¢394H WB »F15e393H —- / 
ese LA »F5e3s4H LBwsF 1503) 
me meRMA TC T-8s14520H MOMENT..Y 1 INEAR WA. »F15.e394H WB »F15e303H = / 
ieee GH wA »F ¥5-3%4H LB »F15-3) 
ORMA ° HtTST LOADS 214) 
2> MeRMAT (T8s20H LIST DISPLACEMENTS 314) 
Pee? 01. Le ls NULDS Z oe ee 
Bee KID COCT) =—1) "20092009 ,2006 
06 GON LINWE = —— cs Lees Se Se Ss ee 
IF (XA-XB) 20083200752008 
00 ONT INVYEW 
RADIUS(1)=XGUEP(IV=YA 
: G@sTOU2004 _. . es 
2008 CONTINUE 
PAD TUS (1) =XGLLD( Tp XA 
GO TO 2004 
2009 CONTINUE = 
Peal S @lebh=SQR FT labXGLieD (1 )-XA)**24+7YGLLDCUI J-YAI**2) 
m 2004 CONTINUE 



































UNIFORM W »F150e394H LA »5F156353H =e 




































































en. OT Won 2 00ms2003,200. " = 
2002 CONTINUE | | a | — = 
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oN AL Re eee UNIO - ———____—___—__— 
STALME(1)=STALD 
GO-TO 203- 
2003 CONTINUE 
4 AEN=SORT( (XB—-XA)**#2+(YB-YA) ¥¥2) 
G INITIAL LOADING IS LINEAR 
ee SA Eel) = S eat DaeeR AD TUS CT AE Ny see Nya Spe 
201 CONTINUE 
—— tl K DAL Dile@alod Viceel Vieni F005 5.350 ee - — 
2005 CONTINUE 
SDp@L? Oe ale 1 sNULDS — = _— —— _ eee: 
DO 210 L=2sNULDS 


Pm ORADTWUS(LI-RADIUS(1)) 20232103210 
mo 202.-CONT I NUE _— ae 5 —— 
NULOM=NULOME(T) 
Ol loeMGLL DC Te —— 
Y6LW=YGLED( I) 
oC LOC AEXLOCAL (TY) 
KINDL=KINDLOCT) 
Pepe LUSRADTUS ( 1) 
STALM=STALME( 1) 
NULOME (CT) =NULOME(L) —_ 
FMoLtLD( 1) =XGLLDI(L) 
een | (ae), = YGieeDeiey A E—E—E—EEEEEEeeee 
POC AL OT V=XLOCAL (L) 
ee UKNDLOCTOEKINDLO(L ) 
Pep lUS(1)=RADIUS(L) 
ome MES IN =STALME(L) 
NULOME (L) =NULOM 
TE < C) 
meek Dtieh=YGLL 
Be OC ALL )=XLOCA 
KINDLO(L)I=KINDL 
| RADIUS(L)=RADIU - 
STALME(L)=STALM 
210 CONTINUE 
meo| SHORTED 
TINUE ic a es 
Oe TO 210252120252102 221035 210232102 32103) »KIOUT 
poze 03 GONZ 1 NUE. P 
DO 2104 T=1»eNULDS 
NKEKITNDOL © a5 
So 10 C7T05 9290622107 %s JUNK 
OAReOnNTLINWE 
Gi PeNeE INTERSECTS A MEMBER 
TIM OL Oe a 
NGAMMA=NUMI JO (MEMBER) 
ND = 
RLEME=SQRT( (HOCOR(NDELTA) ~HOCOR(NGAMMA) )**2+(VECOR(NDELTA) -VECOR ( 
NGAMMA) )**2) : ———— | eg a 
PEPrOX LOGAL 1) -RLEME/2¢6%7 2108 5s 21065 2109 
emmmmrmersie ereT(.T = NU SOC NGAMMA acc yg gr rer 


| NADD=0 
ge ee en hhh 
- 2109 JOINT(1)=NUJO(NDELTA) 
a. hCCMCMCTCUNNAADD =O ee a 
GO TO 2104 
™ 2106 CONTINUE <n _ “a 
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JOTN LOPS NUL OME) eS 
NADD=0 

GOm T Ome? LOM q 
CONTINUE 
MEMBER=NULOME ( L) 
NGAMMA=NUMIJO( MEMBER) 


Bo COND EL TASNUP) JO(MEMBER) 


—E————E————— 


JOINT(1T)=NUJO(NGAMMA) 
JOINT (1+1T)=NUJO(NGAMMA) 
NADD=1 








2 1m, CONMLINUE ai 


NUL=NULDS+NADD 


= tee CK TOUT H4) 27106297 152110 EE EE aaa 


ee OTN TOPS OINT (1-1) ) 211352112 62973 


a 


c 


) 2111 


21% 


ee 15 


ciel 2 


Zig} 


2114 


JOINO=JOINT(1) 
SE ete S ) wal) NO we 
DO 2112 [=2sNUL 











JOINS=JOINT(I) 

; 7 ou) JOTNQ eee 
CONTINUE 
aooe10 360 — 
JOINS=JOINT(1) 
WRITE (7524) JOINS 

Go 2114 1=2sNUu 

me (JOtNT ( [9] JOINT(I-J)) 211552114521 15 
SeiNS=JOINT(I) 

eee FE ee7 24 )= JOIN <2 

CONTINUE 

ie hOesG60  — : — 









































No 2 


Z| 


CONTINUE 
\ Oy SL) 2 isa a 
CONTINUE 


Pee ONE CONGENTRAgWED IwOAD wee 00 EEE 


c 


ee SO led WRATEA 7 +10). 





FACTOR=SQRT ((XA-XB) **2+(YA-YB)**2) 


cee) SeGLORE( (STALD+ENDLD)/Z2ge) — 





GO TO 300 
mmrNwF | sss eal 
MORE THAN ONF CONCENTRATED LOADS 

FACTOR=SQRT(( XA-XGLLD(1) )*¥*2+(YA-YGLLD( 1) )**2) 














Plad=FACTOR*((STALD +STALME(1))/2.-) 
B21) 3 2» NUL DS. 





Zar 








le | 
PACTOR=RADIUS(L)I=-RADIUSC! 
Mer AGTOR*(2-*ABS(STALME(L)) +ABSUSTALMEUI)777(36*( ABS(STALME( LEVIS 





Pee STALME(T)))) ee oe - oe ae 

















P(L)=C#¥((STALME(L)+STALME(1))/2e) 
Be PA | ae G@Re( CSTALME(L.)+STALMEU]) Iyeel Pt) 








Sarr i Nvue 
FACTOR=SQRT((YB-YGLLD(NULDS) )**2+(XB-XGLLD(NULDS) ) **2) 
P(NULDS)=P(NULDS)+FACTOR*((ENDLD +STALME(NULDS) )/2¢) 





eI GON LNWE eee a 


301 


——_ 























THERE ARE NULDS CONCENTRATED LOADS 








Mme (K WP O(K)=1.) 320530) 5 SO 7 ee ce eee 
CONTINUE 

iad. — ~~... si 

IF (1OP=1) 30112301253012 
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31 PF OPATOP +} a : 
GO TO (302530323304) sNLOTY 
4 302 NULOM9=NULOMF (K ) 
P9=P(K) 
a __WRITFE07s121)-NULOM9 sP9 
GOmTO saiz0 
Se 303 NULCMS=NULOME(K) 
P9=P(K) 
ee RoE (7912) NULOMO PS 
GO TO 320 
a 304. NULOMO=NULOME(K) 
P9=P(K) 
es er O71 3) NUL OMO POE ae 
a2 SGONT INUE 
ee JOON) 4OADS “O.Ke= MEMBER LOADS NOW 
IOP=0 
pee) 340 KE) sNULDS “=e 
Mm KINDLO(K)-1) 32133405340 
One \ 
C MEMBER CONCENTRATED LOADS 
emeermen POPs) ) 37115371223712 
alee! TES (7.514) 
meee 12 bOP=)0P +1 = 
Sem O 0322332303247 sNUOTY 
N UMG=fh OM K 
P9=P(K) 
XLOCAS=XLOCAL(K) | 
WRITE (7515) NULOM9sP9»XLOCAY 
oe (Oe TO. 340 sq 
323 NULOM9=NULOME (K) 

































































XLOCA9=XLOCAL (K) 
WRITE (7516) NULOM9sP9+XLOCAQ . ——_ a = 
GO TO 340 
eo NIL OMS=NULOME(K) 
P9=P(K) 
Proc AO=Xi| OCAI (K) $$$ 
WRITE (7817) NULOM9sP9sXLOCAY 
Peon TINUE « — _— 
Cc THIS LOADING ARRANGEDe TRY NEXT 
oe eeorTO 370 
350 GONT INUE 
eel GONTINUE 
a 6 PMBPR t INE t OADS ea ee 
NULD=NULDS+1 
meer ME(NULD)=SENDILD — 
ITOP=0 
bi@ an’ = \ B 
T=Lt+] 
ee] OP = ) 35 1 52 9 3527 " 
Bere WRITE (7514) 
meee 2 1 P= 1OP+) —_— ea 
GO TO (35333543355) sNLOTY 
ON NUE EEE EE 
Me ULTYPE=]) 3537335325353!) 
Sees tL eNULOMO=NULOME(L Jo 
STALM9=STALME(L) 
- ee ed OGAI=XI OCAL CE Dee 
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= YLOCA9=YLOCAL (1) -_ 
WRITE (7918) NULOM9sSTALM9sXLOCA9sYLOCA9 
a GO--7 Omea6 0 c. 
3532 NULOM9=NULOME (L) 
a STALM9=STALME(L) 
STALM8=STALMF(T) 
oe LOC ADS XL OCAR Oe. ee eee 
YLOCA9=YLOCAL(L) 
WRITE (7921) NULOM9 sSTALM9sSTALM8 sXLOCAI»YLOCAD 
ae TO 360 
le 9 4 CON LINUE a 
me CL TYPES FP) Ses 35475 3554) 
mee aes PEL COMO NAL OME Cb a eae 
STALM9=STALME(L) 
POC AQG=XLOGAL (ol) 
YLOCA9=YLOCAL (L) 
See Ee 7 > £9) NULOMOsS TALMO esx! OCADsYLOCAQ 
GO TO 360 


) S PAL Mee 5 PAMME (lm) 
— mime MB=STALME (1)~ 


XLOCAS=XLOCAL (L) 
peel! OCAI=YLOCAL(L) . Sl 
WRITE (7922) NULOM9sSTALM9sSTALM8 sXLOCA9 sYLOCAY 
GO TO 360 
Seo CONTINUE 
eee ee —)) 3551955529355) ee rx orien 
355] NULOM9O=NULOME(L) 
ee Wal MOsSTALME(L) 
XLOCAQ=XLOCAL (L) 
BC AOesGWOGAL (lM) «a EE eee 
WRITE (7320) NULOM9sSTALM9sXLOCA9sYLOCA9Q 
ee TO. 360. 
3552 NULOM9=NULOME (L) 
iia Bee MOeSTALMELI) ———— — 
STALMS8=STALME €1) 
BeAS=X1 OCA — —— 
MmEOCAS=YWOCAL (L) 
WRITE (7923) NULOM9 »sSTALM9sSTALM8sXLOCAIsYLOCAY 
360 CONTINUE 
2 OmeON TINUE —" ge 5. ee — 
DESTROY RADIUS 
wRESTROY STAILME 
BESFROY JOINT 
SweovROY PP _ al | ae 
RETURN 
a UF END a =. 
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me + * O76. _ 
SUBROUTINESAREALO 
DYNAMIC ARRAY  HOCOR sVECOR sNUJOsJORENU sKINDJOsNUPL JO»sNUMI JO »sNUMEMB 
DYNAMIC ARRAY KINDMEsNUJO1sHOCOR1,VECOR]I 
DYNAMIC ARRAY YMEMBsXSTAMEsSTALDMsENDIL DM ~ a 
COMMON HOLENs VELEN» THICK sMANUJO s YOUNG s PRATI OsHOLEN1 sHOLEN2 s VELEN1 
— Od VEL EN2 SEX oWHY SKIND I 9 EX2 so WHY? »PO! OM) oKILILO st TYPE sNUHOU) sNUVE JI » 
2NJsHOCOR(P) sVECOR(P) sNUJO(P) s JORENU(P) sHOLEN3 sVELEN3 sNUHOJ2 s NUHOJ3 
— 3 sNUVEJ2 sNUVEJ3 sNUPLJO(P) »sNUMIJO(P) sNUMEMB(P) »>KINDME(P) »>NUMEsNUVEMB. 
4sNUHCMB sNURIDI sNULEDI sNUPOLOsNUJO1(P) sHOCOR1(P) sVECORI1(P) »sNSIZE> 
Ii SNULDSsENDLDsLDMSTAsLDMENDs KINDJO(P) sK IOUT sNLOTY 
REAL LDMSTAsLDMEND 
ee  PQUITVALENCE. (EXeXA) s (WHY »¥A) s (EXD eXB) of? » YB) seis Tepes POLOM]) ae 
DEFINE YMEMBsNUHOMB+NUVEMB 
= _ DEFINE XSTAME »>NUHOMB+NUVEMB 
DEFINE STALDMsNUHOMB+NUVEMB 
4a _ DEFINE ENDLDM »>NUHOMB+NUVEMB 
Poet YPE-1) 10510320 
® ONT NUE a 
C LOADING UNIFORM OR CHANGING WITH XeDISTRIBUTE AMONG HORIZeMEMBERSe 
n. 1@=0 
DO 100, INDEX=1 »sNUHOMB 
m=} NDEX+NUVEMB | 
NALPHA=NUMIJO(T) 
NA RENE RO ( Ta ES 
IF (XA=HOCOR(NBETA)) 1121003100 
fee CONTINUE : 
IF (XB- —~HOCOR(NALPHA) ) LOUR 1 OOF 





———_ 
















































































Pee CORY NUE nS = ee 
c Se eErLAP POSSIBLE BUT NOT CERTAIN 
= OK NES ‘A 13%)/00 LEE eee 


13 CONTINUE 
i eeieAP=S iidel POSSIBLE 
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PRESET INTEGER *K LOUT + FQUAL 8 —__ ee EE—E—EEE eee 
EXECUTE *OUT! | ee 
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EXDLD 
= = 
Ex2 = 
8OCOR = 
HOCORI = 
ROLEN = 


HOLENT = 


HOLEN2 = 
HOLEN3 = 
JORENU = 
KINDJO = 
KINDME = 
KIOUT = 
Kiibit <= 
LDMSTA = 
LDMEXND = 
LTYPz = 
BANUJO = 
Nd = 


NUJO = 
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APPENDIX E 


A DICTIONARY CF TERMS IN COMNON 





End load magnitude of line loads 
x-coordinate of point supports 

End x-coordinate of line supports 
xecoordinate of joint 

x-coordinate of a non support joint 
Span of plating (parallel to the x-axis) 
Distance along the x-axis from the center of coordinates 
up to the left edge of hele 

Span of hole (parallel to the x-axis) 
NOLEN - HOLENI - HOLEN2 

Joint release number 

Kind of joint 

Kind of member 

A number that specifies necessary output 
Kind of point support 

Start load magnitude for area loads 

End load magnitude for area loads 

A number that specifies kind cf lead 
Maximum number of joints 

Number of joints 

Identification number of joint 


w,..5.P".. a 






















ee eee 
ee eo 
eset bh ene ee Oe 
— ee « 

— ne ee te Oe | 
— << ee) ee ee 
——=———— SF ee ot Gag at ee mi & 
-_— oa Meld ad 
ee ee 
— am . are - 

a mm Dh t 

ae we Ge * 

a 
ee Re 
—— See: 

— ot < Oise es woe 6 
— et ee aes ee 
-——— a, A ome ds 
—- oe 

a & ee - 
el 











ad 














157. 


HUJOL = Identification number of nonesupport joint 

wWULDS = Number of grid loads used to aporeximate a given LINE - 
LOAD 

RUMEMS = Identification mmber of member 

MUMIJO = Identification number of aims joint 


HUME = Sumber of menbers 

HUPLJO = Identification mmber of plus joint 
NULEDI = Number of left diagonal members 
MUPRIDI = Number of right diagonal nenbders 
NUXOMB = Ruamber of horizontal members 

WUVEMNB = Number of vertical menbers 

NOUV EJ 1 

NUVEJ2 

WV AI3 umber of division soints explained in DIFGHO 
NOHO L 

NOHGJ2 

M2053 

KLOTY = Number indicating kind of load 
PRATIC = Poissen's ratio 

POLOML = Point load magnitude 

THICK = Thickness 

VECOR = y-eccord. cf joint 

VECOR1 = y-coord. of non suport joint 

VELEN = Width of plate 

VELEN1L = yecoordinate of lower edge of hole 
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VELEN2 = Width of hole 
VELEN3 = VELEN + VELENL - VELEN2 


WaY 


i 


y~coordinate of a point support 
WiY2 = y-e@nd ccerdinate of a line support 
YOUNG = Young's medulus of elasticity 
KILILO 


= Extra variables not used by prograns 
SUPOLG 





APPENDIX F 





the following instrections are particularly usefrl in using 
the program as it is now, The inmut of course has to be given ace 
cording te FORMATS already included in the progran,. Integers have 
to be right justified in their field, Floating point numbers can 
be anyvhere in their field as long as the point is punched. The 
units to be used ere HIPS and ICES. 

Necessary cards are: 
Gard 1 Forsat (5F 10,3) 
In the first 10 columns put the span of the plate, 
In coluans 11-20 put the width of the plate. 
In columns 21-30 mt the thickness of tne plate. 
In columns 31-40 put Young's sodulus. 
In columns 41-59 out Poisson's ratio. 
Gard 2 Format. (16) 
In columns 1-6, right justified put the desired maximum muabor 
of joints, 





In column 4 put; 


1 if there is no hole 
2 if there is a hole 
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In ease tnere is a hole additional inforsation is necessary in 
this card, Kavely: 
Columns 5-14 Distance along span from the left edge of the plate 
to the left edce of hole. 
Columns 158<4 Distance along width from the lower edge of the 
plate to the lower edge of hole, 
Columns 25-38% Length of hole (in the direetion of the span). 
Columns 35<44 width of hole (in the direction of the width of the 
plste). 

4 Foreat (2 
In coluzns 1-5 risht justified, put the number of point supports. 
In eolums 6-10 right justified, put the mumber cf line supports, 





There must be as many such cards as there are point supports. 
In column & put a: 

1 if support is hinged 

2 if support is fixed 
In cclums 5-14 pot the x-coordinate of particular point support. 
(See Fig, 4,1) (fer location of axes) 
In eclumns 5-14 put the y-coordinate of particular point suppert. 








There must be as many such cards as there are line suppor 
In column 4& put a: 

1 if support is hinged 

2 if support is fixed 


Ait nie nen et! ULL aE ee Ft 
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— 


161, 


In columns $414 mut the xecoordinate of starting point of a line 
support, (In all tne subsequent description the starting point 
of a line is the one to the left if the line is not parallel to 
the y-axis, or the ona downwards if the line is parallel to the 
yoaris) 
in eclumns 15024 put the y-ccordinate of the staring point of the 
line supcort. 
In colunns 25434 put the x-coordinets of the end soint of the 1i- 
ne support. 
In columns 35-4 put the y-ccordinate of the end point of the line 
support. 

6 Fo dy 
In columns 1-4 put the number of point Loads, 
In columns 5-5 put the mumber cf line loads, pight justified. 
In columns 9-12 pat the number of area loads, picht justified. 


In column 5 put as 





1 af lead is a 2 Loree 

2 4f load is an x moment 

3 if load is a y moment 
(The system of coordinates is assumed right-harmied and the loads 
are positive if thoir vector points towards the cositive direction 
of the axes) 
In cclums 6-15 put the x-ccordinate of the point of application 
of the load. 
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In columns 15-25 put the y-csordinate of the point ef applieation 
of the load. 
In columns 26-35 put the magnitude of the load. 


there are as many such cards as there are point loads. 





There are as many such cards as there are line loads. 
In soluen 5 sut a mamber (1,2 or 3) with the same neaning as in 
previous card, 
In eoluane 6-15 put the xecoordinate of the starting point of the 
Line load, 
In colunns 16-25 pit the y-coordinate of the sterting point of 
the lins load. 
In eclusns 26-35 put the magnitude of the load at the starting 
point. 
In column 40 put a: 

0 if line load is uniform 

1 if line iead is linear 
In eoluans 41-50 put the xecoordinate of the end-point of the 
line lead. 
In eoluens 51-60 put the y-coordinate of the and-point of the 
line load, 
In columns 61-70 put the magnitude of the load et the end point. 





There are as many such ecards as there are area loads. 
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In column 5 gut a number (1,2 or 3) with the same meaning as in 
previcus card. 
In columns 6Gel5 cut the x coerdinate of the left lewer corner of 
the sres load, 
in eviumns 16-25 st the y coordinste of the left lower corner 
of tlhe area load. 
In columns 26-35 pat the magnitude of the load at the starting 
point. 
In column: (0 put as 

Q af load is uniforn 

4 if lead is linear along the xedirection 

2 if Load is Linesr slong tie yeiirection 
In eclumns 41-50 put the x-coordinate of the upper right corner. 
In ecluss 51460 mut the y=coomlinate of the upper right corner. 
In eclunns 61~70 pit the magnitude of the load at the end point. 





There are as many such cards as there are kinds of cutput 


desired. 
In column 5 put 2; 
@ 4f all reactions are desired 
3 if all leads are desired at a particular point 
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8h 4£ the loads are desired along a line 

5 if 211 loads are casired 

6 if displanesents are desired at a particular coint 

7 if displacements are desired along a line 

8 if all displacements are desired. 

In case the above mimber is a 2,5 or 8 no other information 
is necessary on this eard. 

In ease the above mmber is a 3 or a 6: 
in ecluens Gel5 wut the x-coordinate of the point where outmt 
is desired. 
In columns 16-25 pit the y-coordinate of the point where outrat 
is desired. 

In case the ahove number is a 4+ or a 7: 
In eclumns 6-15 put the x-coordinate of the starting point of the 
line where ortmt is desired. 
In colurns 16825 put the y-coordinate of the starting point ef the 
line where oautout is desired, 

In columns 26-35 mt the x-eoordinate of the end point of 
the line where outruat is desired. 
in columns 3645 put the y~coordinate of the end point of the line 


where outrmt is desired, 
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‘oLlowing is a preposed user's formulation of a problen, 
the form of COL indieated in Appendix D, is used. 
The statonents have to bo in the order indicated. The under- 
Lived cart of a statenent indicates that this ils the way this state- 
ment is gtcred in the dictionary. Date identifiers can be omitted 
if the deta follow in tho specified ordere 


Any LLST cosmands can appear one after the other. 
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